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ABSTRACT 


The  purpose  of  this  thesis  is  to  examine  the  distribution  of  the 
square  of  the  T  statistic,  under  the  assumption  that  the  data  are 
sampled  from  a  stationary  first  order  autoregressi ve  Gaussian  process. 
For  a  sample  consisting  of  two  observations  the  exact  probability 
density  of  F  =  T2  is  derived,  while  for  large  samples,  asymptotic 
approximations  to  the  probability  density  and  moments  of  F  are 


obtai ned. 
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CHAPTER  I 


INTRODUCTION 


Consider  two  times  series 


X ( t )  =  n(t)  +  e(t) 


and 


Y(  t)  =  v(  t)  +  n(  t) 

for  all  times  t  e  T  ,  where  e(t)  and  rj( t )  are  random  error  terms. 
Let  t^  <  t2  <  ...  <  tn  be  a  finite  set  of  time  points  and  let 
( xi ,y i )  =  (x(tj_)  ,y(tj_) ) ,  (x2,y2)  =  (x(t2)  ,y( t2) ) ,  •  .  •  ,(xn,yn)  = 
(x(tn),y(tn))  be  n  pairs  of  observations  from  the  two  time  series. 
We  wish  to  test  the  null  hypothesis 

H0 :  p(t)  =  v ( t )  ,  t  e  T 

against  the  alternative 

Ha:  ^(t)  *  v(t)  ,  t  e  T 

on  the  basis  of  the  n  differences 

di  =  xi  "  y i  »  1  =  • 

-  1  - 


If  Hq  is  true  then 


1  » 


2. 


e(t.  )  -  Ti(t.  )  =  y. 


1,2 


If,  furthermore,  the  {y.}  are  independent  N( 0 , 1 )  random  variables 
then  the  statistic 


n  d . 

where  d  =  l  —  and 
i  =1  n 


has  a  Student's  t- 


distribution  with  n-1  degrees  of  freedom,  and  the  statistic 


has  an  F-di stribution  with  1  and  n-1  degrees  of  freedom.  However, 
if  the  {y • }  follow  a  stationary  first  order  autoregressi ve  Gaussian 
process 


then  the  distributions  of  T  and  F  are  unknown. 

Chapter  II  surveys  the  literature  dealing  with  the  effects  of 
dependent  data  upon  the  distributions  of  the  T  and  F  statistics. 
Various  types  of  dependence  are  considered. 

In  Chapter  III  we  examine  the  distribution  of  the  F  statistic 
when  the  data  are  sampled  from  a  stationary  first  order  autoregressi ve 


■T  . 


. 

* 


Gaussian  process.  For  a  sample  of  size  n  =  2,  the  exact  probability 
density  of  F  is  shown  to  be 


d>(F) 


/( T-p ) ( l+p) 


tc[1  +  p  +  ( 1 — p ) F]/F 


0  <  F  <  “  . 


It  is  also  shown  that  for  large  n,  the  probability  density  of  F  may 
be  approximated  by 


Knpl/2(1‘pPl) 

/F  ( l-p2p1) 


-(1  +  0(  n-1 ) } 


0  <  F  <  °° 


wi  th 


Pi  = 


_  (l+p2)  +  F( l-p)2/(n-l)  -  V[(1+p2)  +  F( l-p)2/(n-l)]2  -  4p* 


V 


and 


Kn 


(l-P)r(n/2) 


/nii  r[(n-l)/2] 


i-  ( 3  n-1  n  2\  p(n-l)  v(z  n+1  n+2  .  2 

LF  U  •  t  ;  7  ;  p  /  -  s~ir-  •  -T-’  -r-  ;  p  , 


(n-1) (n+1)  J 3  n+3  n+4  .  2] 


n( n+2) 


71  *  ~~T~  ’  ~~T 


p2(n-l) (n+2) (n+3)  c(z  n+5  n+6  .  2^ 

n( n+2)  (n+4)  l\7  *  AT  ’  T"  ’  p  /J 


-1 


where  F(a,p;y;z)  is  the  hypergeometric  function.  Finally,  large 
sample  approximations  for  the  moments  of  F  are  obtained. 


CHAPTER  II 


LITERATURE  REVIEW 

If  _X  =  (Xj_ ,X2» . . .  ,Xp) '  is  a  vector  of  n  independent  identically 

distributed  normal  random  variables,  each  with  mean  \i  and  variance 
2 

a  ,  and  i f 


then  the  statistic 

T  =  x  -  P 
S//n 

has  a  Student's  t-di stribution  with  n-1  degrees  of  freedom,  and 

c  _  t2  _  n(X-n)2 

r  -  i  -  k- 

has  an  F-di stribution  with  1  and  n-1  degrees  of  freedom. 

Independence  of  the  {X^.}  is  a  major  assumption  used  in 
determining  the  distributions  of  the  T  and  F  statistics.  Thus  it  is 
of  considerable  importance  to  examine  the  sensitivity  of  these 
statistics  to  conditions  of  dependence.  In  this  section  we  review  the 
relevant  literature. 

Let  {en}  be  a  sequence  of  statistics  such  that 

9n  -  E(9n> 

rmj 


5. 


is  asymptotically  distributed  as  a  standard  normal  random  variable,  and 
consider  testing  the  null  hypothesis  HQ:  E(e  )  =  0  .  The  critical 
point  corresponding  to  a  right-tailed  test  of  size  a  is 
determined  by  the  relation 

1  -  $(K  )  =  a  , 
a 

where  $  is  the  standard  normal  distribution  function.  If  we  denote 
by  D  the  asymptotic  variance  of  9n  in  the  case  of  dependent 
observations  from  a  stationary  process,  its  value  in  the  case  of 
independent  identically  distributed  observations  by  V,  and  the  ratio 
D/V  by  x  ,  then,  as  shown  by  Gastwirth  and  Rubin  (1971),  the  critical 
region  determined  under  the  assumption  of  independent  observations, 
namely 


/TTe~T 

has  approximate  probability 

a*  =  1  -  ®(K  x'1/Z)  . 
a 

Following  these  authors,  we  define  a*  to  be  the  asymptotic  level  of  the 
right-tailed  test  based  on  en  when  the  observations  are  dependent. 

Scheffe  (1959)  considered  the  model  where  X.  =  (X^ ,X2, . . . ,Xn) 1  has 
a  multivariate  normal  distribution  with  mean  vector  jj,  =  (p,p,...,^ a)1 


. 


6. 


and  covariance  matrix  £  ,  where 


I  = 


1 

P 


p 

1 


p 

1 


(0) 


(0) 


1 

P 


P 

1 


To  ensure  that  l  is  positive  definite,  he  imposed  the  necessary  and 
sufficient  restriction 


-^2 [cos 


n 

TTTTrn 


-l 

} 


<  p  <  < 2  cos 


t 


'  it 


-1 

) 


Under  these  conditions,  he  showed  that  T  is  asymptotically 
N ( 0 ,  l+2p)  ,  and  that  the  asymptotic  level  of  the  right-tailed  T-test 
of  Hq :  [i  =  0  is  given  by 


a*  =  1  -  m  (l+2p)‘1/2) 
a 


Let  _X  =  (X1,X2, . . .  ,Xn) '  have  a  multivariate  normal  distribution 
with  mean  vector  0_ =  (0,0,...,0)‘  and  covariance  matrix  l  .  Let 
X  =  and  let  denote  the  element  in  the  ith  row  and  jth 

column  of  \  •  Hotelling  (1961)  demonstrated  that  for  large  t,  the 


n"/2Ul1/2 


(I  l 


Wl 


stati Stic 


- 


approximates  the  ratio 


where  I 


he  found 


while  for 


with  p  > 


P(T>tl X-1  =  \  ) 

p<T>tir1  =  i) 


is  the  n  x  n  identity  matrix. 


For  the  covariance  matrix 


n-1 

p 

n-2 

P 


•  •  • 

n-1  n-2  n-3 

P  P  P 


Rn  to  be 


Rn  ■ 


(1-p) 


1/2 


(l-p)n[l+2p(l-p)-1n'1]n/;! 


P  P  •  *  *  P 

1  p  .  •  •  P 

P  1  «  •  •  P 


•  • 


P  P  P 


-(n-1)"1  ,  he  showed  that 


5 


R. 


-ri  + 


(n-l)/2 


8. 


In  this  latter  case  Ali  (1973)  proved  that  the  exact  probability 
density  function  of  T  is 


9(t) 


r(£) 


r(^)U(n-l))1/2(l+  y^)1/2 


TrTTT  +  H 


1 

TT 


W2 


He  also  showed  that  the  tails  of  the  t-di stribution  for  p  >  0  are 
thicker  than  those  of  the  standard  t-di stribution  (p=0)  ,  while  for 
p  <  0  the  tails  are  thinner  than  for  the  standard  case. 

Gastwirth  and  Rubin  (1971)  examined  the  behaviour  of  the  asymptotic 
levels  of  the  sign,  Wilcoxon,  and  t-tests  for  various  dependent 
processes.  They  showed: 

i)  If  X^,X2»...,Xn  are  from  a  completely  regular  Gaussian 
process  whose  autocorrelations  satisfy  l  ~  and 
p^  >  0  for  all  k,  then  the  level  of  the  sign  test  is  less 
sensitive  to  the  dependence  than  the  level  of  the  Wilcoxon 
test  which  is  less  sensitive  that  of  the  t-test.  In  each  case 
the  level  exceeds  the  corresponding  value  for  independent 
observations. 

ii)  If  X1,X2,...,Xn  are  from  a  first-order  autoregressive 

Gaussian  process  with  negative  autocorrel ation  coefficient 
p  ,  then  the  level  of  the  sign  test  changes  less  with  the 
dependence  than  the  level  of  the  Wilcoxon  test  which  changes 
less  than  that  of  the  t-test. 


9. 


Albers  (1978)  studied  the  model  in  which  the  random  vector 
X_  =  ( X-j^ , X2 , .  - . , Xp )  has  a  multivariate  normal  distribution  with 
vector  \x  =  (  (j.,p, . . .  ,|i) 1  ,  and  the  {Xn*  }  are  m-dependent;  that  i 
covariance  matrix  l  has  elements 

a..  =  Cov  ( ,Xj  )  =  a2P|i_j.  |  ,  1  <  I  i  -  j  I  < 

=  0  ,  I  i  -  j  I  >  m  , 

for  i , j  =  l,...,n,  where  m  is  a  positive  integer  and  , 
k  =  l,...,m  are  constants  such  that  £  is  positive  definite  and 

m 

1  +  25;  p  >  0  . 

k=l  K 

Under  the  given  model,  he  showed  that  Tf  is  asymptotically 

N(0,  1+  2  l  p  (1-  t)) 
k=l 

and  that  the  asymptotic  level  of  the  right-tailed  T-test  of 
Hq:  p=0  is 


m 


a*  =  1  -  ®(Ka{l  +  2  |  pk(l-  ^)}-1/2)  +  0(1) 

k  1 


mean 
,  the 


This  same  paper  also  dealt  with  the  situation  in  which  the  {Xn-}  are 
from  a  stationary  autoregressi ve  process  of  order  m;  that  is, 


5 


i  =  m  +  1 , . . . ,  n 


10 


m 

I  3.  ( X .  .  -|i)  =  2  • 
k=0  K  1_K  1 


2 

where  zm+1,...,zn  are  independent  N(0,y  )  random  variables,  a0  =  1 
and  ai,...,am  are  constants  such  that  all  roots  of  the  equation 


lie  inside  the  unit  circle.  Under  these  conditions  he  showed  that  T 
is  asymptotically 


m  /  m  2 

N(0,  l  a  Pk/(  l  ak)  +  0(£))  , 

•k=0  K  V  k=0  K  n 

where  is  the  autocorrel ation  coefficient  at  lag  k,  for 

k  =  0,1, ... ,m. 

For  each  model,  Albers  proposed  a  robust  modification  of  the  t- 

test.  He  demonstrated  that  under  independence  both  of  the  proposed 

2 

tests  would  require,  asymptotically,  mk^  additionally  observations  to 
obtain  the  power  of  the  t-test. 


CHAPTER  III 


THE  ASYMPTOTIC  DISTRIBUTION  OF  THE  SQUARE  OF  THE  T  STATISTIC 
FOR  A  FIRST  ORDER  AUTOREGRESSIVE  GAUSSIAN  PROCESS 


3.1  First  Order  Autoregressi ve  Processes 

Consider  a  finite  subset  of  time  points  t^  <  t2  <  ...  <  tk  and 
let  F[X(tj) , X ( 1 2 ) >• • • , X ( t k ) U  denote  the  joint  distribution  function  of 
the  random  variables  X(t^),  X( t 2 ) » • • • ,X(tk ) .  A  time  series  or 
stochastic  process  is  said  to  be  strictly  stationary  if 

F[X(t1),X(t2),...,X(tk)]  =  F(x(t1-x),X(t2-T),...,X(tk-x)] 

for  all  nonempty  finite  subsets  t-j.  <  t2  <  . ..  <  tk  and  all  x  ;  that 
is,  if  the  joint  distribution  is  independent  of  the  time  origin. 

A  process  is  defined  to  be  weakly  stationary  if  for  all  t,  t*, 
and  x 

1)  E[X(t)]  =  E[X(t-x)], 

ii)  V[X(t )]  =  V[X(t-x)]  ,  and 

iii)  Cov[X(t )  ,X(t*)]  =  Cov[X(t-x) ,X(t*-x)]  . 

That  is,  for  a  weakly  stationary  process  the  mean  and  variance  are 
independent  of  time,  and  the  covariance  is  invariant  under  a  change  of 
time  origin. 

Consider  the  first  order  autoregressive  process 
(3.1.1)  X^.  =  pX-£_i  +  e -j-  9  t  =  0,  +1,  +2,... 


11 


. 


12. 


where  the  {et}  are  uncorrelated  random  variables,  each  with  mean  0  and 
o 

variance  a  .  Following  Fuller  (1976),  it  may  be  seen  by  repeated 
application  of  (3.1.1)  that 


xt  =  P2xt-2  +  Pet-1+  et 

=  P3xt-3  +  P2et-2  +  Pet-1  +  et 

!  .  k-l 

=  P^xt-k  +  l  P1et  • 
i=0 

Assume  that  | p|  <  1.  Then 


almost  surely,  and  thus 


(3.1.2) 


E(Xt) 


E(  I  p\)  =  0 

i=0 


and 


Cov(Xt,Xt+h)  =  E[(Xt-E(Xt))(Xt+h-E(Xt+h))] 


=  E(XtXt+h)  =  E 


L  Yj  =  0 


l 

P  e 


t-i 


PJe 


t+h-j 


2  r  i  i+h 
=  cr  l  p  p 

i  =0 


* 


>  K  qq 


13. 


almost  surely.  Thus,  for  |p|  <  1  the  first  order  autoregressive 
process  (3.1.1)  is  weakly  stationary.  Suppose  in  addition  that  the 
{e^}  are  uncorrelated  N(0,a  )  random  variables,  so  that  any  finite 
subset  of  the  (X^}  has  a  multivariate  normal  distribution.  Then, 
since  a  multivariate  normal  distribution  is  determined  by  its  first  two 
moments,  the  modified  process  is  strictly  stationary. 


. 


14. 


3.2  The  Joint  Distribution  of  X-|  ,X^> 5 . . . 

Let  X^,X2,...,Xn  be  sample  observations  from  the  stationary  first 
order  autoregressi ve  process 

(3.2.1)  Xj.  =  +  e^  ,  t  =  0,+l  ,+2  ,... 

where  | p|  <  1  and  the  {et}  are  independent  N(0,1)  random 
variables.  From  equation  (3.1.2)  we  have 

h 

E(Xt)  =  0  and  Cov(Xt,Xt+h)  =  . 

1-p 

Thus  _X  =  ( Xj ,X2 , . . . ,Xn ) '  has  a  multivariate  normal  distribution  with 
probability  density 


f(XJ  =  (2it)"n/2|£r1/2exp 


(X'  l'h) 


where  l  is  the  nxn  covariance  matrix 


1 

P 

i — i 

i 

c= 

CL 

C\J 

T 

2 

•  •  •  r\ 

1-P 

1-P 

1-p 

CL 

1 

n-2 

P 

1-P2 

i-P2 

•  •  •  r\ 

1-P 

• 

• 

• 

• 

• 

• 

• 

• 

• 

n-1 

P 

n-2 

P 

1 

i-P2 

i-P2 

•  •  •  r\ 

1-P 

15. 


The  following  method  of  evaluating  \l\  is  due  to  Patton  (1961) 


m  ■ 


1 

CL 

1-P2 

1-P2 

CL 

1 

1-P2 

• 

1-P2 

• 

• 

n-1 

P 

n-2 

£ 

i  2  i  2 
1-p  1-p 


n-1 

p _ 

1-p2 
n-2 
£ _ 

i-P2 


i-p' 


Multiplying  the  second  column  by  p  and  subtracting  from  the  first 


m  - 


1 

0 

• 

p 

n-1 

P 

1 

P 

n-2 

P 

r\  •  •  • 

1-P 

1 

i-p2 

n-2 

P 

i-p2 

P 

r\  •  •  • 

1-PZ 

1 

i  3 

1-P 

n-3 

P 

r\  •  •  • 

1-P 

• 

i-p2 

• 

i  2 

1-P 

• 

r\  •  •  • 

1-P 

• 

i-p2 

• 

• 

0 

n-2 

P 

• 

1 

n-2 

P 

n-3 

P 

• 

1 

r\  •  •  • 

1-P 

i-p2 

n 

i-p2 

r\  •  •  • 

1-P 

i-p2 

n-1 


Continuing  in  this  fashion  we  find  that 


m  ■ 


i-p' 


. 


16. 


We  wish  to  determine  an  expression  for  .  Following 

Routl edge's  (1972)  derivation  of  ,  we  define  the  nxn  matrix  J 

by 


J  = 


1  -p 

1  -p 

1  -p 


(0) 


(0) 


1  -p 
1 


Then  l  =  j£  is  the  nxn  lower  triangular  matrix 


I  = 


.  p 

2 

1 

1 

(0) 

P 

p 

• 

• 

• 

• 

• 

• 

• 

• 

• 

n-2 

n-3 

n-4 

...  1 

P 

p 

P 

n-1 

n-2 

n-3 

p 

P 

P 

P 

1-P2 

i  2 

1-p 

1-p2 

•  •  •  r\ 

1-P 

To  obtain  the  nxn  identity  matrix  I  from  l  we  perform  the  elementary 


column  operations  -pC^  + 


-pC3  +  C2  -pCn  +  Cn_1  , 


(1-p  jC  .  Applying  these  same  linear  transformations  to  l_  we  see 


that 


(0) 


1 

-p  1 

-p  1 


(0) 


-p 


1 

-p 


Thus 


/V 

r1  ■  r1! 

=  1  -p 

-p  1+p2  -p  (0) 

-p  1+p2  -p 

•  •  • 

•  •  *9 

(0)  -p  1+p 

-p 


and 

x'r1^  =  X1  +  +  ( 1+P2) (x|+. . 


-  2p(X1X,+X,X,  +  ...+X„  ,Xn  ,+Xn  ,X 
12  23  n-2  n-1  n-1 
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3.3  The  Exact  Probability  Density  of  F  for  a  Sample  of  Size  n  =  2 
In  this  section  we  determine  the  exact  probability  density  of  the 
F  statistic  for  the  case  in  which  the  sample  consists  of  2 
observations . 

Let  and  ^  be  sample  observations  from  the  stationary  first 
order  autoregressi ve  process  (3.2.1).  Then  and  X2  have  a  joint 

bivariate  normal  distribution  with  mean  vector  0_  =  (0,0)'  ,  covariance 
matrix 


1 

p 


p 

1 


and  probability  density 


f(x1#x2) 


a7 

2k 


exp 


'|(X^-2pX1X2+X^) 


In  addition  we  have 


where 


-  Xx+X2  2  (X1-X)2+(X2-X)2  (Xx-X2)2 

X  =  — ^ —  and  S  = - 2 — ^  —  2 

Thus 


F 


<  W*  _  u2 

(X1-X2)2  V2 


where  we  have  introduced  the  transformati on 
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(3.3.1) 


U  =  *1  +  X2  ,  V  =  X-l  -  X2 


The  inverse  transformation  of  (3.3.1)  is 


X,  = 


U+V 


1  2 


y  -  u’v 
X2  '  ~ 


with  -oo  <  u  <  ®  and  -<*>  <  V  <  oo  .  The  joint  probability  density  of 
U  and  V  is  given  by 

g(u,v)  =  fCX^UjV)  ,X2(u,v)]  I J  I  , 

where 


J  = 


d(xx,x2) 

1/2 

1/2 

5(u,v) 

- 

1/2 

-1/2 

-1 

7 


Thus 


9<u,v)  =  exp^ 


-  2p 


_ (-1  fu^(  1  — p )  +V ^ (  1+p) 


Tt~  exp\  7 


■]} 


which  is  the  probability  density  of  a  bivariate  normal  random  vector 
Y_  =  ( U , V ) '  with  mean  0_  =  (0,0)'  and  covariance  matrix  with 


l nverse 


-I 


kf-  0 


0 


1+p 


By  inspection  we  obtain 


0 
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1-p 


2 

T+p 


5 


and  therefore  U  and  V  are  independent  normal  random  variables  with 


and  V 


It  follows  that 


^pLu  and 

are  independent  standard  normal  random  variables  and 

U2  and  U^iv2 


are  independent  chi-square  random  variables,  each  with  1  degree  of 
freedom.  Thus  the  statistic 


F* 


u-p!  uz  / U+p)_  v2 


.  (i-p)  F 

'  TT+pT 


has  an  F-di stribution  with  =  1  and  v2  =  1  degrees  of  freedom. 

Let  4>  and  p*  denote  the  probability  densities  of  F  and  F*, 


respectively.  Then 


. 


g.  bfti  -  fttiw  £6(1 


p*(F*)  = 


v  i /2  n  1 

r[(v1+v2)/2](v1/v2)  1  F*vi 

r(v1/2)r(v9/2)  (Vl+v9)/z 

1  ‘  (i+v1f*/v2)  1  ^ 

r( l )  F*"1/2 

=  r(  l/2)r( l'/'Zy  TI+PT 

=  1 
%[  1+F*  )/FF 


and 


<D(F)  =  d>*(F*(F)) 


dF* 

dF 


=  _ 1 _  (1-p) 

it[  l+(  1-p  )F/  ( 1+p )  ]/( 1-p  )F  /  ( 1+p )  ^  p  ^ 

Thus  the  exact  probability  density  of  the  F-statistic  is 


*( F)  =  /(1~P)(1+P) 


■jiE  l+p+(  1-p )  F]/F 


3.4  The  Cramer-Geary  Inversion  Formula 

Daniels  (1956)  outlined  a  proof  of  the  Cramer-Geary  inversion 
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formula.  The  following  derivation  is  based  upon  that  presented  by 
Patton  (1961). 


Consider  a  statistic  of  the  form 


where  Cq  is  almost  surely  positive.  Let  f(cQ,c^)  denote  the  joint 
probability  density  of  c0  and  Cj_ .  We  wish  to  obtain  the  distribution 
of  r  from  f(cQ,c^). 

The  transformation 


c0  -  c0  >  C1  "  rc0 

has  Jacobian  * 


d(Cg,C^) 

m 

1  0 

a(c0,r)  " 

r  c0 

so  that  the  joint  probability  density  of  Cq  and  r  is 


c0f(c0,rc0) • 

Integration  over  Cq  yields  the  marginal  density 

oo 

(3.4.1)  h(r)  =  /  c0f(c0,rcQ)dc0  , 

which  is  the  probability  density  of  r. 


Let  M(Tq,T1)  denote  the  joint  moment  generating  function  of  c0 
and  cj_  .  Then 


WTq,^)  =  E  \e 


T0C0+T1C11 


00  CO 


,  ,  TOcO+Tlclf,  X  ,  , 

=  /  /  e  f  (cA,c1  )dcAdc 


-CO  0 


o’  r  0  i 


' 
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By  the  Fourier  inversion  formula,  it  follows  that 


ffCQ.C^  = 


( 2ni ) 


/  /  MdQ.T^e 


-(T0c0+tici) 


dTodTi 


The  integration  is  along  the  imaginary  axes  in  the  Tq  and  planes 
from  -i®  to  +i®  ,  or  along  any  allowable  deformation  of  these  paths; 

that  is,  along  any  paths  from  -  i®  to  ^  +  i®  such  that  no 
singularities  of  M(Tq,T^)  lie  on  the  new  paths  of  integration,  or 
between  them  and  the  imaginary  axes.  Substitution  of  rcg  for 
gives 


1  -Cndo+rT,) 

(3.4.2)  f(c0,rc0)  =  -  //  M(Tq ,T1  )e  dTQdT1  . 

( 2ni.) 

Under  the  linear  transformation 

u  =  To  +  rTi  ,  Ti  =  Ti  , 

with  Jacobian 


d ( Tq  *T 1 ) 

l 

-r 

a(u  ,Lr1)"" 

0 

1 

equation  (3.4.2)  becomes 


f<VrC0>  =77^//  M(u-rTi,T1)e  0  dudT1  . 

(Z-rtl  ) 

The  integration  of  u  is  along  a  path  in  the  u-plane  correspond!" ng  to 
that  of  Tq  in  the  Tg-plane.  Thus 
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I 

0 


c  u 

f ( c q , rc q ) e  0  dcQ 


=  / 


0 


(  1  _cnu  )  cnu 

— \  2  //  M(u-rT1,T1)e  dudl^fe  dc 

(  (2ni)  1 ) 


0 


Ml 


'  l  ‘cnu 

/  M(u-rT1  ,T1  )e  u  du 


cQu 
e  dco[dT1 


By  the  Fourier  inversion  formula. 


Mfu-rT^Tj) 


/ 

0 


•  1  -cnu 

^  /  M(u-rT1,T1)e  du 


C0U  , 

e  dcQ  , 


so  that 


(3.4.3) 


c~u 


/  f(cQ,rc0)e  0  dcQ  =^-/  M(u-rT1,T1)dT1  . 


When  differentiation  under  the  integral  sign  with  respect  to  u  is 
permissible,  we  have 


c^u 


/  f(c0,rc0)c0e  dcQ  f  5u 


0 


M(u-rT1,T1) 


dTL  . 


Letting  u  =  0  ,  we  obtain,  by  equation  (3.4.1) 


h(r)  =  "ZiW  5u  CM(  u-rTi  ,Ti )] 


dT, 


u=0 


Often  it  is  desirable  to  transform  to  another  variable  z 
by  T-l  =  T^(z,u),  where  T].(z,0)  maps  the  T^-plane  onto  some  region  of 
the  z-plane.  In  this  case 


aT. (z,u) 

M( u-rT  ,T  )dT  =  M[u-rT  (z,u)  ,T  (z,u)]  — — -  dz 

111  i  i  oz 


If 
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which,  for  notational  ease,  we  express  as 


M(u-rT1,T1) 


dz  . 


From  equation  (3.4.3)  it  follows  that 


00  cou 

/  f(cn,rcn)e  dc 

0 


1 

0  2n;i 


/  M( u-rT x , T x ) 


dz  . 


The  integration  is  along  the  transformed  path  in  the  z-plane.  Once 
again  differentiating  with  respect  to  u  under  the  integral  sign  and 
letting  u  =  0,  we  obtain 


(3.4.4) 


h(r)  /  It 


M(u-rT,  ,T.) 


5T- 


1*1  dz 


dz 


u=0 


This  equation  is  known  as  the  Cramer-Geary  inversion  formula  and  will  be 
used  in  the  work  to  follow. 


•i*. 
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3.5  The  Joint  Moment  Generating  Function  of  Cq  and 

Let  Xj_,X2» . .  •  ,Xn  be  sample  observations  from  the  stationary  first 
order  autoregressi ve  Gaussian  process  (3.2.1)  and  let 


where* 


X 


and  S' 


I(XrX)2  IX2  -  (£x.)2/n 
— rwl  =  rTT 


Then 


where 

n(IX. )2 

n£X2-(£x.)2 

Let  g  and  h  denote  the  probability  densities  of  F  and  r, 
respectively.  Then  g  and  h  are  related  by 


(3.5.1) 


and 


(3.5.2) 


g(F)  =  h(r) 


dr 

-  n  h 

nF  1 

dF 

TFIT  h 

(n-l)  

h(r)  =  g(f) 


dF 

W 


(n-l)  On-l)  r 
n  L  n 


*  All  summations  in  section  3.5  are  from  i  =  1  to  i  =  n  unless  otherwise 
i ndicated. 
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We  may  express  r  as  the  ratio 


where 


and 


We  shall  use  the  Cramer-Geary  inversion  technique  to  obtain  an 
asymptotic  approximation  of  h(r),  from  which  we  shall  approximate 
g(F).  In  this  section  we  derive  an  expression  for  M(Tq,T^),  the  joint 
moment  generating  function  of  Cq  and  c^. 

By  definition. 


(3.5.3)  MCTq.T^  =  E(e 


n 


where 


x'rh  =  X1  +  xn  +  ( 1+p2)  (x|+. . -+xn_i 
-  2p(X1X,+X,X,  +  ... 


12  2  3 


. 


. 


Si  nee 


1  w  r-1, 


Vo +  Vi  -  4 x n 


=  (X^+X^)[(n-l)T0+nT1-  |]  +  (X2+. .  .+X^_1)[( n-l)TQ+nT1  -^(1+p2)] 


+  (XlX2+X2X3+"-+Xn-2Xn-l+Xn-lXn)[-2T0+2nTl+P] 


+  (X.Xo+X^.  +  .-.+X^  +X,X.+X,Xc+...+X,X  +...+X„  ,X  )[-2Tn+2nT.]  , 
1  o  1  4  1  n  d  4  d  b  2  n  n -d  n  u  1 


we  have  that 


(3.5.4) 


Toco  +  Tici  -  ^x'r1*'  =  *^x'bx) 


where  B  is  the  nxn  matrix 


B  = 


f 

b 


b 

a 

b 


b 

a  b 


(c) 


(c)  b  a  b 

b  f 

with 

a  =  1  +  p2  -  2(n  —  1)Tq  -  2n T , 
b  =  2Tq  -  2nT1  -  p  , 
c  =  2Tq  -  2nT1  ,  and 
f  =  1  -  2(n-l)T0  -  2nT1  . 


Combining  equations  (3.5.3)  and  (3.5.4),  we  obtain 


,  .  % 
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dX 


which,  by  the  theory  of  the  multivariate  normal  distribution,  reduces  to 


(3.5.5) 


From  equation  (1.23)  of  Appendix  I  we  have 


n-1 


<j[u+(  f-a)]2[(  2e+d)  2-c( ( u-v)-n( 2e+d) )] 


I B I 


(u-v)(2e+d)  1 
+c( 2e+d)[u2-( f-a) 2] j 

- - - 7i[v+(f-a)]2[(2e+d)  2+c(  (u-v)+n(  2e+d) )  ] 

( u-v) ( 2e+d) L  ' 

+c( 2e+d)[v2-( f-a) 2]| 


+ 


(2e+d)c 

where  u  and  v  are  roots  of  the  equation 


-  dx  +  e^  =  0 


wi  th 


and 


e  =  b  -  c  =  -p  . 
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This  expression  for  |B|  is  symmetric  with  respect  to  u 


Also,  the  roots  of 

x^  -  dx  +  e^  =  0 


are 


+  v42-4e' 


so  that 


(3.5.7) 


,2  ,  .2  ,  2, 
d  -  d  -4e  ) 
uv  =  - -i - - 


=  e 


and  v  . 


3.6  An  Asymptotic  Approximation  for  | B | 


Setti ng 
(3.6.1) 


u 


we  obtain,  from  equations  (3.5.6)  and 
-  dz“l  +  e^  =  0 


and 


v  =  e2z  =  P2z  . 

Let  |i  =  Tq  +  rT^  .  Then  from  equations  (3.6.1)  and  (3.5.6) 

z~^  +  e2z  =  d 


or 


(3.6.2)  z_1  +  p2z  =  1  +  p2  -  2nTQ  =  1  +  p2  -  2n^ 
When  p.  =  0  equation  (3.6.2)  becomes 

(3.6.3)  z'1  +  p2z  =  1  +  p2  +  2nrT1 


and 


or 


P  z  - 


=  - 
2nr|_  z 


(3.6.4) 


(l-z)fl 
2nr  Lz 


The  mapping  (3.6.3)  consists  of  the  mappings 
1  2 

A:  z_i  +  p  z  =  w 
and 

2 

B:  w  =  1  +  p  +  2nrT^ . 


(3.5.7) , 


we  have 


2nrT^  . 


p2( 1-z) 


. 
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A  is  a  simple  modification  of  the  Joukowski  transformation 

z“l  +  z  =  w  , 


which  is  discussed  in  Appendix  II.  By  arguments  similar  to  those 
presented  in  Appendix  II,  it  may  be  shown  that  A  maps  the  w-plane,  cut 
along  the  real  axis  from  -2 1 p I  to  2|p|,  inclusive,  onto  the  interior 

of  the  circle  |z|  =  -Jv  . 

Ipl 


Under  the  second  mapping,  if  w  =  -2 1 p  I  then 


T  _  (1+2  Ip  I +p  )  _ 

'1  "  2nr 


(1+lpl) 

2nr 


and  if  w  =  2 | p |  then 


Ti  ■  - 


( 1-2 1 p |+p^)  _  ( 1- I p  I 


2nr 


2nr 


Thus,  B  maps  the  T^-plane,  cut  along  the  real  axis  from  - 


to  -  ,  onto  the  cut  w-plane.  For  Ipl  <1  and  all  positive 


n  and  r,  the  endpoints  of  the  cut  in  the  T^-plane  are  negative. 

Of  particular  interest  is  the  transformed  path  of  integration,  that 
is,  the  path  in  the  z-plane  into  which  the  imaginary  axis  in  the  T^- 
plane  is  transformed.  Let  z  =  x  +  iy.  Then  equation  (3.6.4)  becomes 


_  ( 1-x-iy )/x-iy 
1  -  ^%2+y2 


l-x)-iy 


x 


L  x  +y 


2 

P 


“27T 

x  +y  . 


5 


so  that 
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Re(Tl>  ■  TW 


(l-x)l^^- -  p2U  y 


2  , 


2  2 

x  +y 


~2  2 
x  +y  J 


The  imaginary  axis  in  the  T^-plane  has  equation  Re(T^)  =  0  and  thus 
on  the  transformed  path  of  integration 


0  .  X  2  X2  2  y2 

0 - * — -  p o — o  +  Xp  - 


~2  1 

x  +y 


2  2 

x  +y 


2  2 
x  +y 


It  follows  that 


=  x|  -p-~y  +  p^)  -  (1+p^)  • 
^x  +y 


(3.6.5) 


x^  +  y^  = 


1+p  ( 1 -x ) 


and 


-V 


-  X  =  + 


1+p  (1-x) 


x( 1-X ) -p2 ( l-x)x2 

1+P2(1-X) 


x ( 1-x ) ( 1-p  x) 

1+P2(1-X) 


Since  y  is  real,  we  have 


x ( l-x)( 1-p  x)  >  Q 
1+p2 ( 1-x) 


But  I  pi  <1  and  |z|  <  y^y-  ,  and  so  we  must  have  0  <  x  <  1  .  We 


also  see  that  for  any  given  x,  lyl  is  a  decreasing  function  of  p‘ 
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which  is  maximized  when  p  =  0  .  In  this  case  |y|  =  /x( 1-x)  ,  which 
has  maximum  value  1/2  when  x  =  1/2.  Thus  0  <  |y|  <  1/2  on  the 
transformed  path  of  integration.  From  equation  (3.6.5)  it  follows  that 

I  z  |  =  /x2+y2  =  ^  / £ -  ,  0  <  x  <  1 

\lV(l-x) 

which  is  an  increasing  function  of  x  and  has  maximum  value  1  when 

x  =  1.  Thus,  I z I  <  1  on  the  path  of  integration.  Additionally,  we 

see  from  equation  (3.6.4)  that  as  T^  ->  0,  z  +  1  and  as  T^  -►  +  i®  , 

z  +  0.  Therefore  as  T^  traverses  the  imaginary  axis  from  -  i®  to 

0,  z  moves  clockwise  along  a  path  in  quadrant  I  from  (0,0)  to  (1,0), 

and  as  T^  goes  along  the  imaginary  axis  from  0  to  +  i®  ,  z 

travels  clockwise  in  quadrant  IV  from  (1,0)  to  (0,0). 

1  2 

Substituting  e  =  -p  ,  u  =  z  ,  and  v  =  p  z  into  the  expression 
(3.5.6)  for  I  B_| ,  we  see  that  the  second  term  (that  is,  the  term  with 

v11"1)  involves  a  factor  p2(  n-1 )  zn-l  and  the  term  (that  is,  the 

term  with  (-e)n_1)  involves  a  factor  pn~'1'  .  Since  Ip  I  <1  and 
| z |  <  1  *,  it  follows  that  both  terms  will  make  contributions  which  are 
exponentially  small  for  large  n.  Furthermore,  in  the  case  p  =  0  they 
make  no  contribution  at  all.  Thus  we  may  neglect  these  terms  and 
approximate  |B|  by 


*  In  fact  | z I  <  1  on  the  transformed  path  of  integration  except  at  the 
point  z  =  (1,0). 
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(3.6.6)  | B |  -  - H - J[u+(f-a)]2[(2e+d)2-c((u-v)-n(2e+d))] 

( u-v) ( 2e+d) L  ' 

+  c(  2e+d)[u2-(  f-a) 

=  u.n..-i— [u+(  f-a )  ][  ( 2e+d)  2-c  ( ( u-v )  -n(  2e+d ) )  ] 

( u-vH  2e+dp ' 

+  c(  2e+d) [u-f+a] J-  . 


In  the  special  case  (p=0)  ,  (3.6.6)  is  an  exact  expression  for  |B|. 
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3.7  The  Special  Case  (p=Q) 

Later  we  shall  use  the  Cramer-Geary  inversion  formula  to  obtain  the 
asymptotic  distribution  of  the  r  statistic.  To  illustrate  the  method, 
we  derive  the  exact  distribution  of  r  for  the  special  case  (p=0). 

Let  p  =  0.  From  equations  (3.5.4),  (3.5.6),  and  (3.6.6)  we  have 

I B I  =  u°  Lu+f-a)  j[u+( f-a) ][( 2e+d^)-c( (u-v)-n( 2e+d) ) ] 

( u-v ) ( 2e+d  j 

+  c( 2e+d)[u-f+a]j 


i  2  2 

where  u  =  z  ,  f  -  a  =  -p  =  0,v  =  pz  =  0,e  =  -p  =  0,  and 

i  2  -1 

d  =  z-i  +  p  z  =  z  .  Thus 


IB!  ■  --TV' 


TV 


) 


{ 


z"1[z”2-c(z”''-nz’1)]  +  ctz”1 


Hz"1)) 


01* 

(3.7.1)  IBI  =  -i^l  +  ncz}  . 


Solving  for  T-^  in  equation  (3.6.2)  we  obtain 
(3.7.2)  Jl  =  ^iFCz_1-l+2np.]  . 


From  (3.5.4)  we  have 

c  =  2T0  -  2nJl  =  2p  -  2(n+r)Tx 

so  that 
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c  .  2,  - 


=  2(i  - 


2(n+r)|i  (n+r)  (1-z) 


nr 


-2n|i  (n+r)  (1-z) 


nr 


It  follows  from  equation  (3.7.1)  that 


IBI  =4{i  .2Dj£.i"i!±(1.2)}  , 


and  so,  by  equation  (3.5.5), 


(1-P2)172 


M(T0’T1>  =^T72 


IBI 


|B| 


1/2 


=  zn/2{l  -  2p2^  -  (n+r)  (i_z)| 


-1/2 


Differentiating  both  sides  of  equation  (3.7.2)  wih  respect  to  z 


we  obtain 


dT 


1  -1 


az  2^7 


so  that 


MfTo.Tp 


aTj  _z(n-4)/2 


dz 


2nr 


2n2^tz  (n+r)  _1//2 


(1 


->> 


and 


- 


Thus 

(3.7.3) 


where  z  =  x 
i ntegrati on 


2n 


_nZ(n“2)/2ji  (n+r) 

— ^2— i1  -  — 


(1-z)} 


-3/2 


-nz^n"2^/2((n+r) r  r 


2r 


t 


r  (n+r) 


-  1  +  z 


} 


-3/2 


nz(n~2)/2  L  .  r 
2/F(n+r)3/2(  ln+r) 


} 


-3/2 


J,.jU 

2|z  Tn+rTf 


-3/2 


2/nr( 1+r/n) 


h  ( r ) 


- - - 377  !  1 

4-n:i/nF(  1+r/n) J/  r 


dz 


P=0 

(n-2)/2J 


\z  “  Tn+rl 


-3/ 


4n;i/nr(  1+r/n)3^2  r 


/  z(n-2)/2(z-p)-3/2dz 


+  iy  ,  0  <  p  =  n/( n+r)  <  1  ,  and  r  is  the  path  of 


r:  lyl  =  /x( 1-x )  , 


0  <  x  <  1 
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which  is  a  circle  of  radius  1/2  centred  at  (1/2,0). 

For  p  =  0  ,  we  know  that  the  F  statistic  has  an  F  distribution 
with  =  1  and  vg  =  n  -  1  degrees  of  freedom.  It  follows  that  F 
has  probability  density 


g(F)  = 


r[ ( v^+V2)/2] ( v^/ Vg) 

r(v1/2)r(v272)" 


v-|  /2 


v i7  2— ! 


( 1+v^F/ v  £ ) 


( v^+v2)/2 


=  r(n/2)[l/(n-l)]1/2  F'1/2 

r(  1/2  )r[  ( n-1  )/2 )  J  j-  /  ( n-1 ) 

=  r(n/2) _ 1  _ 

/(n'-rr^F  r[(n-l)/2]  [  1+F/ ( n- 1 ) ] n/ 2  ’ 


From  equation  (3.5.2)  we  have 


(3.7.4)  h(r)  =  g[— 


-1)' 


n  J 


.  (n-1) 


r(n/2) 


n  /(n-1  )tc( n-1  )r/n  r[(n-l)/2]  {1  +(n-l)r/[n(n-l) ]} 


WZ 


r( n/2) 


TTZ  ’ 


/n nr  r[(n-l)/2]  (1+r/n) 

which,  to  verify  our  method,  we  must  obtain  from  equation  (3.7.3). 
Consider  the  path  of  integration  r  which  we  may  express  as 


follows : 


r.  z  =  1  +  1  J9 


Starting  from  z  =  0  ,  as  z  describes  the  circle  in  the  clockwise 
direction  9  ranges  from  -n  ,  through  an  angle  of  2tc  ,  to  +n 
Letti ng 


(z-(3)  =  R0e 


i(9+2kn) 


baft 


•; 


40. 


we  have 

( z-p ) “3/2  .  r-3/2  e-i ( 30/2+3kn) 
which  yields  two  solutions 

(z-p)"2//2  =  r;3/2  e"30"1'^2  (principle  determination) 

and 

(z-p)~2/2  =  R”3^2  e"( 3e/2+3n;)i  (secondary  determination). 

0 

-3/2 

To  make  (z-p)  single-valued  we  take  it  to  be  the  principle 

determination  and  cut  the  z-plane  from  -®  to  p  so  that  the  secondary 
determination  is  unable  to  cross  the  cut.  This  procedure  also  makes 
z(n-2)/2  single-valued,  because  in  cutting  the  z-plane  from  -®  to 
p  ,  it  is  automatically  cut  from  -«  to  0. 

Consider  the  closed  path  P  consisting  of  the  paths  r‘  ,  L^,  y  , 
and  L2  defined  by 

1  1  i  fi 

r':  Z  =  Y  +  ?  e  ,  n  -  £  >  0  >-n  +  £ 

L^:  z  =  p  +  Te1^"1x+£\  r(e)  >  i  >  6 

i  0 

y:  z  =  p  +  5e  ,  -it  +  e<9<it-e 

l_2:  z  =  p  +  xeiU  e)  ,  6  <  t  <  r(e)  , 

where  e  is  an  arbitrarily  small  angle,  r(e)  is  the  distance  from 
( p , 0 )  to  the  point  of  intersection  of  and  r  ,  and  6  is 
arbitrarily  chosen,  such  that  0  <  6  <  mi n( p ,1-p)  .  P  is  depicted  in 
Figure  1. 


' 


41. 


Figure  1:  The  path  of  integration  P 


' 


Since  z^  n~2)/2(  z-p)  is  single-valued  and  analytic  on  and 

within  the  closed  path  P  we  have 

(n-2)/2 


/  M  M  } 

Li  y  lJ  (z-p) 


37T 


dz  =  0  , 


so  that 


(n-2)/2 


/ 

r'  (z-p) 


W 


dz  = 


(n-2)/2 


'4  Y  L2  (z-p) 


372 


dz 


On  l_i  ,  dz  =  e1^"7r+£^  dx  ;  on  l_2,  dz  =  e^’Mx  ;  and,  on 
y  ,  dz  =  die1'0  de  .  Thus 


( n-2 ) / 2 


/ 

r'  (z-p) 


3/7 


dz  =  - 


c  it  +  \  (n-2)/2 

6  ( p+xe  l(^+e))  J(-*+e)  , 

— - 77^5 - e  dx. 

r(e)  /  ^i(-n+e) 


( xe 


) 


^  / 


n~£  (p+6ei9)^n~2^2  _ie 


-u+e  (6ei0)3/Z 


6ie  de 


+  / 

5 


r(e) 


(p+xe 


•/  t  ( n-2 ) / 2 
l  C  tc— e )  \ 


) 


i ( n-e ) 


/  i  ( n-e )  J/Z 
( xe  ) 


dx 


! 


■-(I 


5  (fs+zei{-n+c)) 


(n-2)/2 


r(e)  t37y^K77T 
(n-2)/2 


dx 


it-e 


(p+6e19) 


j.  -s  r  vp^uc  / 


.  fr(e)  (p+TeiU-£)) 

+  l  - 3/2 


x  e 


( n-2 ) / 2 
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dx> . 


v 


‘‘v 


' 
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Letting  e  -►  0  and  I  =  / 


( n -2 ) / 2 


r  (z-p) 


tti 


dz  we  obtain 


-■{/ 


^  (p+xe  l71) 


( n -2 ) / 2 


%  ,  „ .  c. _  i  e 


(n-2)/2 


p  x3/V1lt/z 


dx  +  i/  l£i6-e-  ) 

J  TO 


-it  (Se1^ 


de 


o  ,  (n-2)/2 

X  (P+^e  11 )  d  \ 

+  J  - Trr-rm  QT| 


6  xW^" 


/  ^  /  a  n-2)/2 


71  (p+6e10) 


(n-2)/2 


t/5  -n  (ei9)1/2 


de 


P  /Q  x(n-2)/2  \ 

-  '!  “t! -  *} 


or 

(3.7.5) 


I  =  i 


t2I>  •  i  ‘4 


where 


h  - 1 


p  (3-.)(n-2)/2  ,  T  .  r%  (3+5ei6)(n'2)/2 

372  d%  d  *2  ^  “W? 

6  x  —ic  e 


Consider 


-  p(n-2)/2/P^  _  .yn-2)/2x-3/2  dt 


.  Then  x  =  pu  and  dx  =  p  du  ,  giving 

I.  -  p^"2)/2  /  (l-u)(n-2)/2(3u)-3/2P  du 

6/p 


de  . 


Let  u  =  x/p 


<  a 


=  p<n-3'/2  /  (l-u)(n-2)/2  u'3/2 

6/p 


du 


Applying  the  general  binomial  formula,  we  see  that 


y  (  nk  r(n/2) .  k 

k^0  (  1  rck+i) rcn/2-k ) 


converges  uniformly  to  ( 1-u) ^ n~2)/2  for  q  <  u  <  1  . 


Thus 


T  -  0(n-3)/2  v  ,  ,,k  r( n/2) 

1  '  P  k=0  WTTTlnT^  6/ 


/  uk'3/2du 


.  „(n-3)/2  "  ,  ,,k  r ( n/2 ) 

P  L  '  A '  TTTT+TTrTn 


r„k-l/2 


k=0 


r(k+l)r(n/2-k ) 


u 


LTW?) 


] 


6/p 


and  so 

(3.7.6)  21,  -  2e("-3)/2r(n/2)  r(k+l()r(n/2-k)(k-l/2) 


(n-2)/2 

-  r( n/2 )  l 

JE  k=0 


(-l)k(6/p)k 

r(k+l)r('n7'2'-V)'(k-l/2T  ' 


Now  consider 


i2  .  P(n-2)/2/TC 

-n 


ie 


d n-2 ) / 2 


-i  9/2 


de  . 


1  f) 

Since  |6e  /pi  <  1  ,  we  have,  by  the  general  binomial  formula,  that 


v  r( n/2)  ( 6  Je\k 

r(k+i)r(n/2-k)^  ) 


converges  to 


(n-2)/2 


uniformly  in  9  .  Therefore 


■ 
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I2  =  p(n'2)/2r(n/2) 
»  p(n‘2)/2r(n/2) 
=  p(n-2)/2r(n/2) 


Y  (6/p)k  i e( k-1/2 ) 

k;0  rik+i)r(n/2-k)  1 


;  (6/P)k  rei9(k-1/2)  1 

k^0  f(V+l)r(n/2-k)  L  i  ( k'-T/2T \ 


7X 


“IX 


00 


I 

k=0 


( 6/p) k  2 

r(k+l)r(n/2-k)  (k-1/2) 


si n[n(k-l/2) ]  . 


\c  + 1 

Noting  that  sin[^(k-l/2)]  =  (-1)  ,  we  obtain 


(3.7.7) 


2p(n-2)/2 

/6 


r(  n/2) 


"  (-l)k(6/p)k 

kiQ  r ( k + 1 )  r ( n  / 2 - k )  ( V-'172T  ' 


Combining  equations  (3.7.5),  (3.7.6)  and  (3.7.7)  we  have 


x  ,  2i3(n-2)/2 
/6 


r(n/2) 


T  (-l)k(5/p)k 

Lk|0  r(k+l)r'(n/2-k'RF-T7?T 


v  (-l)k(6/p)k 

k|0  r(k+l)r(n/2-k)(k-l/2) 


CD 

+  2i  ^n_3^/2r(n/2)  l 


k=0 


(-IT 

r(k+l)r(n/2-k)(k-i/2) 


(n-3)/2 


=  21 


n+r 


K 


n 


where  Kp  =  r(n/2)  ^  r(k+l)r(n/2-k)(k-l/¥T  ’  From  equation  (3'7'3)  we 


obtain 


46. 


h(r)  = 


1 


4ui/nr( 1+r/n) 


377 


2u/nr( 1+r/n) 


1  /n+r 

377  \~ 


r(n-3)/2 


K 


n 


and  thus 


(3.7.8)  h(r)  =  " - i -iL  . 

/urn  (l+r/n)n/  2/u 

Comparing  equations  (3.7.4)  and  (3.7.8),  we  see  that  it  remains  to  show 
that 


Kn  _  -  r(n/2) 

ITT  rlTn--T)72T  ' 

Renormalization  yields  this  result  since  h(r)  is  a  probability  density 


function. 


47. 


3.8  An  Asymptotic  Approximation  for  h(r). 

In  this  section  we  obtain  an  asymptotic  approximation  for  h(r)  for 
the  general  case  0  <  |p|  <  1.  From  equations  (3.5.4),  (3.5.6)  and 
(3.6.6)  we  have 

IBI  «  u  1(u+f-a),{(u+f-a)[(2e+d)2-c((u-v)-n(2e+d))]  +  c( 2e+d) (u-f+a) } 

( u-v) ( 2e+d) z 

1  2  2 

where  u  =  z_i  ,  f  -  a  =  -p  ,  v  =  p  z  ,  e  =  -p  ,  and 
d  =  z'1  +  P2z  .  It  follows  that 

u  +  f  -  a  =  z"1  -  p2  =  (l-p2z)z"1  , 

u  -  f  +  a  =  z"1  +  p2  =  (l+p2z)z-1  , 

u  -  v  =  z_1  -  p2z  =  ( 1-pz) ( 1+pz ) z” 1  , 

and 

2e  +  d  =  -2p  +  z-1  +  p2z  =  (l-pz)2z_1  . 

Thus 

un-1(u+f-a)  =  z'|n~ll(l-p2z)z-1  =  (l-p2z) 

( u-v)  ( 2e+d)Z  (l-pz)(  1+Pz)z"1[(  l-pzl^z'1]2  zn'3(l-Pz)5(l+Pz) 

and 

|  B |  *  ( l-p2z)  l-p2z)  |-(  1-pz)4  _  c  A  1-pz)  ( 1+pz)  _  n  ( l-pz)2\ j 

zn"3(l-pz)5(l+pz)^  z  7  V  z  1  ' 

+  c  ^  1-PZ7  (  1+pz)2) 

Z  z  ( 


- 

T 
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(1-P  z) 

zn'1(l-pz)4(l+pz) 


c(l-pz)(l+p2z)\ 


3 

|(l-p2z)[-^-'-P-z-L  .  c( l+pz-n(  1-pz) ) ] 


By  equation  (3.5.5),  we  have 


M(Tn,T 


0’ 1 1)  = 


(1-P2) 1/2 


1 B  | 


T7T 


’  (1-p2)1/2[^T 


1 ke!?l _ ]{d-P2z)[(1~pz)3 


Z  '  (l-pz)  (1+pz)  <■ 


-  c( 1+pz-n+npz) ] 


c( 1-pz) ( 1+P2z)|] 


-1/2 


From  Equation  (3.6.2)  we  obtain 


(3.8.1)  T,  =  4br  Cz'1  +  P2z  -  1  -  P2  +  2np]  =  ^  t*1'21!,1'13  z)  +  2np] 


1  2nr 


so  that 


aT) 

5Z 


2nr 


( -z-2+p2 )  =  - 


(1-P2Z2) 

- 5 - 

2nrz 


and 


aTl 

L5Z  J 


=0 


49. 


(1-P2z2)(l-P2>1/2' 


2nrz' 


(1-P  z) 


zn_1( l-pz>4( l+pz) 


{U- 


p2z) 


(l-pz) 


c( 1+pz-n+npz) 


C(l-pz)(l+p2z)l  "3/2( -1/2) 


_  2 

■  j1"9  Z\ - (-(l-p2z) ( 1+pz-n+npz)  +  ( l-pz) ( l+p2z)) 

Lzn“i(l-pz)4(l+pz)  J 


dc 

5p 


p=0 


Combining  equations  (3.5.4)  and  (3.8.1)  we  find  that 


c  -  2Tq  -  2nT^ 


=  2\i  -  2(n+r)T. 


=  2pi  - 


2(n+r) 

2nr 


(l-zMl-p^z) 


(n+r)  ( 1-z) ( l-p2z) 
nr  z 


2n|i 

r 


Therefore 


(n+r)  (l-z)(l-p2z) 
nr  z 


and 


50. 


(3.8.2) 


dT, 


Ldz  J 


1+ 


=0 


d-pVjd-p2)1^ 
“  2 - 


4nrz 


(1-Pcz) 

L2n'1(l-pZ)4(l+pZ) 


{( l-p2z)[(-1-~.ez.! —  +  iHili  U-zHl-p  Z)  ( 1+pZ.n+npZ)  ] 


-  (n^rl  (1-Z)(1-P2Z),  (1.pz)(1+p2z)| 


-3/2 


(1-p  z) 

_zn_1( l-pz)4( 1+pz) 


-  ( 1-P2z)( l+pz-n+npz)+( 1-pz) ( l+p2z) 


=  -  (1-%  z>(1-p  ^  4  (1~p-  z--  [-( l-p2z)  ( 1+pz-n+npz)  +  ( 1-pz)  ( l+p2z)  ]Q  , 
2r  z+i( 1-pz) 4( 1+pz) 


where 


(1-P2z) 

^n“1(l-pz)4(l+pz) 


(1-P2z) 


U-pz)3  + 


(n+r)  ( 1-z) ( l-p2z) 
nr  z 


( 1+pz-n+npz) 


(n+r)  ( 1-z) ( l-p2z) 
nr  z 


( 1-pz) ( l+p2z) 


Consi der 


51. 


r-n,(1-p2f)2  {(1-PZ>3 

Lzn(l-pz)4(l+pz)  ' 

-^-(l-z)(l-pz)(l+P2z) 


+  — ^-■■■(l-z)(l-p2z)(  1+pz-n-npz) 


(1-P2z)2 


_zn(l-pz)  (l+pz)nr^ 


|nr( 1-pz ) 


n(  l-p2z) ( 1-pz)  -  (l-pz)(l+p2z) 


(n+r) ( 1-z) [( l-p2z) ( 1+pz) 


2  2 

Observing  that  (1-p  z)(l+pz)  -  ( 1-pz ) ( 1+p  z)  =  2pz(l-p)  ,  we  obtain 


Q  = 


(1-P2z)2 


zn(l-pz)4(l+pz)nr 


|nr( 1-pz ) 


+  ( n+r)  ( l-z)[2pz(  1-p)  -  n(  l-p2z)  ( 1-pz)  ]  J> 


-3/2 


Neglecting  2pz(l-p)  in  comparison  with  n( 1-p  z) ( 1-pz)  we  incur  an 
error  which  is  relatively  0(n“l),  and  thus 


(1-P2z)2 


zn(l-pz)^(l+pz)nr 

(1-p2z)2 

Lzn(l-pz)3(l+pz)r  _ 


|nr(  1-pz)' 
1  -3/2 


-  (n+r) ( l-z)n(l-p2z)(l-pz)| 


-3/2 


where 


52. 


L  =  r(l-pz)2  -  (n+r)(l-z)(l-P2z) 


=  [r  -  (n+r)]  +  z[-2rP  +  (n+r)(l+p2)]  +  z2P2[r  -  (n+r)] 
=  -n  +  z[n( 1+p2 )  +  r(l-p)2]  -  np2z2  . 


Letting  and  p2  denote  the  roots  of  L  =  0,  we  may  express  L 

2 

L  =  -  np  (z-p1)(z-p2)  ,  so  that 


(3.8.3) 


-nP2(l-p2z)2(z-p1)(z-p2)  " 
-  zn(  1-pz)  l+pz)r 


-3/2 


where  0  <  |p|  <  1  , 


Pi  = 


n( 1+p2)  +  r( 1-p2) 


-  A n(l+p2)  +  r(l-p)2]2  - 

2nP2 


4n2p2 


and 


q  _  n( 1+p2)  +  r( 1-p) 2  +  /^n ( 1+p2)  +  r(l-p)2]2  -  4n2p2 

Let  us  digress  briefly  to  examine  the  roots  and  p2  .  The 
expression  under  the  radical  sign, 

[n{ 1+p2)  +  r< 1-p)2]2  -  4n2p2 

=  n2 1 1+p2) 2  +  2nr( 1+p2) ( 1-p) 2  +  r2( 1-p)4  -  4n2p2 

=  n2(l-p2)2  +  2nr( 1+p2) ( 1-p)2  +  r2(l-p)4  , 


as 


53. 


is  real  and  positive.  Thus  both  p1  and  p2  are  real  and  p2  >  0  . 
Furthermore, 

Q  o  _  Cn(l+p2)  +  r(l-p)2]2  -  { [n ( 1+p2)  +  r(l-p)2]2  -  4n2p2} 

'2  ^77 - 


2  2 

_  4nV  _  1 

"  ? 
4n  p  p 


and  so,  >  0  .  Also,  since 


n2(l-p2)2  +  2nr( 1+p2) ( 1-p)2  +  r2(l-p)4 


=  n2( 1-p2) 2  +  2nr( 1-p2) ( 1-p)2  +  r2( 1-p)4  +  4nrp2(l-p)2 


>  n2( 1-p2) 2  +  2nr( 1-p2) ( 1-p) 2  +  r2( 1-p)4 


=  [n( 1-p2)  +  r(l-p)2]2 


we  have  that 


„  _  n( 1+p2)  +  r(l-p)2  -  /n2(l-p2)2+2nr(l+p2)(l-p)2+r2(l-p)4 

Pi - 2 - 

1  2nP^ 


„  n( 1+p2 )  +  r(l-p)2  -  /[n(l-p2)  +  r(l-p)2]2 

<  - - 

2np^ 


_  2np  _  , 

- ? - 1  * 

2nP^ 
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Thus  0  <  p  <  1  and  p2  =  1/ ( P* 2p1)  >  1/p2  >  1 


From  equations  (3.8.2)  and  (3.8.3),  we  obtain 


—  M( T  T  ) 

dp.  U  O’V 


'dL 

dF 


=0 


-  -  (1-/i)(1-pY/2(1'p2z)  [~( l-p^z)  ( l+pz-n+npz)  +  (l-pz)(l+p2z)] 
2r  zn+i( 1-pzj  ( 1+pz) 

-np2(  1-P2z)  2(z-p1)(z-p2)  _,~3//2 


zn(l-pz)3( l+pz)r 


] 


=  (12PnZl)(1-p241/2(1~p2z)  [(1-p2zH1+pz)  -  n(l-p2z)(l-pz)  -  ( 1-Pz) ( l+p2z) ] 
2r2zn  i(l-pz)4(l+pz) 


'--nP2(l-p2z)2(z-p1)(z-p2)'l_3/2 


zn(  l-pz)3( l+pz)r  J 


Once  again,  we  oberve  that  by  neglecting  the  terms 

2  2  2 

(1-p  ) ( 1+pz )  -  ( 1-pz) ( 1+p  z)  in  comparison  with  n(l-p  z)(l-pz)  we 

introduce  an  error  which  is  relative  0(n“M.  Thus 


p: 


=0 


M  2  2Ul  2xl/2M  2  v  9 

(  T  Ll  P  4  U~P-  '-[-n(l-p2z)(l-p)] 

2r  z+  ( 1-pz ) 4 ( 1+pz) 


“-z^(l-pz)3(l+pz)r 

jip2(l-p2z)2(z-p1)(z-p2). 


3/2 


which  upon  simplification  becomes 
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(3‘8‘4)  k*{W 


■&T. 


dz 


i(1,p2)l/2z(n-2)/2(1,pz)5/2(1+pz)3/2 


lp=0  2p/nr(l-Ptz)(z-p1)^'(z-p2) 


It  follows  that 


(3.8.5)  h(r)  = 


1  r  *— 

7nT  Jrie  dU 


r  5Ti" 

M(M,-rT1,T1)  — 


dz 


P=0 


1  ,  i(l-p2)1/2Z(n-2)/2(l-pz)5/2(l+Pz)3/2 

7*r J 


r*  2p3/nr  (l-p2z)(z-p1)3/2(z-p2)::1/2 


.  U-P2)1/2 


/ 


f(z) 


4np3/nr  r*  (z-p1)3/2 


dz 


dz 


where 


f(z) 


:<n-2)/2(l-pZ)5/2(l,pz)3/2 

(l-p2z)(z-p2)3/2 


z  =  x  +  iy,  I z |  <  1  ,  and  r*  is  the  path  of  integration 


r*  :  lyl  = 


\ 


x( 1-x) ( 1-p  x) 


l+p7(l-x) 


From  the  analysis  of  the  path  of  integration  presented  in  section  3.6, 
we  know:  (i)  r*  is  a  closed  path  which  is  symmetric  about  y  =  0; 
(ii )  0  <  x  <  1;  (iii)  0  <  lyl  <  1/2;  (iv)  as  Jl  travels  along  the 
imaginary  axis  from  -i°°  to  0,  z  moves  in  a  clockwise  direction  in 


. 
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quadrant  I  from  (0,0)  to  (1,0)  ;  and  (v)  as  goes  along  the 
imaginary  axis  from  0  to  +i°°  ,  z  proceeds  clockwise  in  quadrant  IV 
form  (1,0)  to  (0,0).  The  path  r*  is  shown  in  Figure  2. 

Consider  the  integrand 


f(z) 


3/2 


2(n-2)/2(i.pZ)5/2(i+p2)3/2 

( l-P2z)(z-p1)3/2(z-p2)3/2 


of  h(r).  There  are  branch  points  at  z  =  p^  ,  p^  ,  1/p ,  -1/p  ,  and, 
if  n  is  odd,  at  z  =  0.  Letting 


(z-^) 


R(0 


^i  ( 9+2k7i:) 


we  have 


(z-Pl)-3/2  -  [R(9)Pl)]-3/2e-i(30/2+3k^  , 
which  gives  two  solutions 

(z-p^-3^2  =  [(R(0,p1)]"3//2e_1  30//2  (principle  determination) 


and 


(z-p1) 


-3/2 


=  [R(9,P1)]'3/2e‘(3e/2+3,l)l 


(secondary  determination) 
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z-pl ane 


0 


Figure  2: 


The  path  of  integration  r* 
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-3/2 

To  make  (z-p^)  single-valued  we  take  it  to  be  the  principle 
determination  and  cut  the  z-plane  from  -®  to  p^  so  that  we  are 
unable  to  cross  the  cut  and  reach  the  second  determination. 

Similarly,  to  ensure  that  z ( n~ 2 ) / 2  n-s  single-valued,  we  take  the 
principle  determination 

z(n-2)/2=  [R(9,0)]<"-2)/2eie("-2)/2  . 


Note  that  in  cutting  the  z-plane  from  -«  to  p^  ,  it  is  automatically 


cut  from  -«  to  0,  as  required.  Also,  to  make  (1+pz) 
valued  we  take  the  principle  determination 


3/2 


si ngl e- 


(l+pz)3/2  =  [R[e  ,-1/p  )]3/2e1' 30/2  ; 


the  z-plane  is  already  cut  from  -®  to  -1/p  . 


Finally,  we  wish  to  be  able  to  cross  the  x-axis  between  p^  and 

1.  Noting  that  1/p  <  1/p2  <  P2  ,  we  cut  the  z-plane  from  1/p  to  =□ 

5/2  -3/2 

and  take  the  secondary  determinations  of  (1-pz)  '  and  (z-p2) 


Taki ng 


(l-pz)5/2  -  [R(e,l/p)]5/2el(50/2+5ll) 


and 


-3/2  _  r„/„  „  , -i— 3/2„— i  ( 39/2+3n) 


(z-p2)‘J/i  =  [R( 9  ,p2)]  e 
we  make  ( 1-pz ) 5/2  and  ( z-p- ) _3/2  single-valued,  and  avoid  the 
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singularity  of  the  integrand  at  z  =  1/p2  . 

Consider  the  closed  path  P*  consisting  of  the  paths  r*' ,  L*  ,  y*  , 
and  L£  def i ned  by 


r*  :  r*  ,  n  -  e  >  9  >  -n  +  s 


Lf :  z  =  $i  +  Tel("1I+£)  ,  r{z)  >  1  >  5 


Y*:  z  =  ^  +  Se 


ie 


-ix+£<9<7i;-e; 


L|:  z  =  (3^  +  t1*  ^ £  ^  ,  6  <  t  <  r(e) 


where  e  is  an  arbitrarily  small  angle,  r(e)  is  the  distance  from 

(p^,0)  to  the  point  at  which  L*  and  r*  intersect,  and  6  is 

arbitrarily  chosen  such  that  0  <  6  <  min( p.  ,1-p. )  .  P*  is  illustrated 

-3/2 

in  Figure  3.  On  and  within  P*,  f(z)(z-p^)  is  single-valued  and 

analytic.  Thus 


v  *  '  •  J  one 


. 
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Figure  3:  The  path  of  integration  P* 


61. 


On  LJ  ,  dz  =  e1^‘7l+£^dx  so  that 


f(z) 


/ 

LJ  (z-^) 


Tfl 


6  f(p1+xei(-7t+£))ei(-1T+e) 

dZ  "  1(e)  (xe1^')3'2  dT 

_  g-i  (  -it+e)/2  J6  ^i+'ce1<  11  e)> 

r(e) 


17T 


dx 


As  e  -►  0  ,  r(  e)  -►  ^  ,  and 


f(z)  dz  +  eU/2  f6f_!h^£2dx. 


s  - 

Lf  (z-px) 


37T 


P 


l/T 


1 


Since  e1*7^  =  cos  %/ 2  +  i  sin  tc/2  =  i  and 

e"1'71  =  cos ( -ic)  +  i  sin(-Tc)  =  -1  ,  we  see  that  in  the  limit  as  e  +  0 


(3.8.7) 


f(z) 


4  (z-Pi> 


3/2 


dz  =  i  / 


6  f(Pj-x) 


3/2 


dx 


Similarly,  along  L|,  dz  =  e1'^71  dx  and 


f(z) 


4  (z-pi) 


TT 


,  x  r/ .  .  i ( n-e) \  i ( u-e) 
r(e)  f(p,+xe  )e 

iZ°l  (xe1*™*)^  ^ 

.-1 (u-e)/2  fr(e)  f(Pl+Te  U  '' 


dx 


Observing  that  e  lix^  =  cos ( — -rc/2 )  +  i  sin(— n/2)  -  -i  and 
e1*71  =  cos  ( it)  +  1  sin  it  =  -1  >  we  obtain,  in  the  limit  as  e  ->  0  , 


f(z) 


4  (z-p^ 


3/2 


Pi  f  (  P i -x ) 

dz  =  -  i  /  - 372”  dx 


(3.8.8) 


. 


* 


- 


62. 


i  f) 

On  y*  >  dz  =  i^e  de  and  so 


f/7\  Ti-e  f (p1+6e1'9)6e1'9i 

/  1  3/2  dz  =  /  —-■■■■«,  ?/y—  de  > 

Y*  (z-^)372  -u+e  (6ei0)2/2 


which,  in  the  limit  as  e  -►  0  becomes 


H  =  —  /  f(p1+5ei9)e-i9/2d9  . 

/6  -71 


Consider  the  circle 


C: 


z  =  ^  +  6e19 


2 

with  centre  ((3^,0)  and  radius  5  .  Since  ^  +  5  <  1  <  1/p  <  ,  C 

2 

does  not  contain  the  points  (1/p  ,0)  or  (p2,0)  .  Thus  f(z)  is 
analytic  on  and  within  C,  and  we  may  expand  f(p^+6e  )  in  a  Taylor 

series  about  .  It  follows  that 


n 


H  = 


—  I 

Jb  -h 


k=0 


f (k) (pL) 

in 


ie,k 


e-l9/2  de 


=  4  I 


/6  k=0 


,k 

-m - 6 


(6e  ) 

i9(k-l/2)  K 


TIF-172T 


—  Tl 


1  °° 

1 

/6  k=0 


f (k) (g  ) 

r  rJic(k-l/2)  rt-in(k-l/2)  i 

k !  ( k'-T/TT  L  J 


Since  e1  ^(k-l/Z)  _  e-inke-iit/2  _  (_1)k(_1*)  and 
e-in(k-l/2)  _  g-i-rikgin/2  _  we  have 


■ 


H  = 


2i  “  f(k)(?!)5k  k+1 

^  J0kT(k-l/Z)  (-l) 
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Noting  that  for  k  >  1 


/  x1^2  dx 
0 


,k-l/2 


(FT/'2) 


1 

k-1/2 


0 


we  obtain 


H  = 


2i_ 

/s 


»  f(k)(P1)6k  ...  1  k 

2f 0,)  +  I  - pi - (-l)k  1  /  xk'3/2dx 

1  k=l  K-  0 


(k) 


1  oo  1 

—  f(pJ  +—  /  -x"3/2  j;  -T1 — —  (-x6)kdx  . 

/6  /6  0  k=l 


But  the  Taylor  series  expansion  of  f(^-x6)  about  is 

»f(k)(M  k 

f(Pi-x6)  =  f(pj  +  l  - n -  (-x6)  dx  , 

i  i  k=1  k. 


and  so,  in  the  limit  as  e  -*  0 


I 


f  ( z ) 


Y*  (Z-0X) 


1372 


dz  = 


—  fop 


/6 


1 


.  21  /V3/2[f(p,-x6)  -  f(p.  )]dx 
/SO  1  1 


(3.8.9) 


Observing  that  as  e  -*•  0  ,  r*'  -►  r*  ,  and  combining  equations 


(3.8.6)  through  (3.8.9)  we  obtain 


(3.8.10) 


/ 


r*  (z-px) 


f(z)  Pi  ^(pi-^)  4.: 

riZ)37?  dz  =  21  /  1  —  dx  -  —  f($, ) 


6  X 


1 


/6 


1 


+—  /  x-3/2[f(3l-x6)  -  f(p.)]dx 

/6  0 


Consi der 


H  =  21  / 

5 


372" 


dx 


Letting  y  =  x/p^  ,  we  have  x  =  p^y  ,  dx  =  p^dy  , 


and 


H,  =  —  J  y'3/2f(p,-Piy)  dy 
1  6/p1  1  1 


2i 


,-3/2 


=  —  I  y_0/  ^[f(p1-p1y)-f(p1  )]dy  + 

6/Pj. 


2if(p.)  1  ,/2 

- L.  /  y'3/2dy 

6/^l 


Now 


/  y‘3/2[f(Pi-Piy)  -  f(Pi)]dy 

6/Pi 


1 


=  /  y"3/2[i:(p1-P1y)  -  f(f3i)]dy  -  h2 
o 


6/pl  -3/2 


H2  =  /  y'  '  Cf(P1-P1y)  -  f ( ) ]dy  . 


0 


where 


.  o 


Substituting  x  =  p^/6  we  see  that  y  =  6x/px  ,  dy  =  (6/p^dx  ,  and 


1 


6x>  -3/2, 


=  /  (^)_"C[f(Pr6x)  -  f(p.)] 

0  P1  1  1 


✓FT  1  o /p 

=  -zr  I  *  y  L f(p,-6x)  -  ftp.)]  dx  . 
✓6  0  1  1 


Furthermore , 


1 

I 

6/Pi 


-1/2 


1 

[l  -  /Pll 

=  -2 

6/Px 

L  /6  J 

=  2 


^1 

/6 


-  2  , 


and  thus 


Hi  ■ 


2i_ 


J  y_3/2[f (p1-31y)  -  f(p:)]dy 


✓FT  i  -3/p 

— -  /  X  J/^[f(p1-6x)  -  f(p1 )]dx 
✓6  0 


2if(p.  ) 

4-  1 

f2^i  1 

J  /Pi 

_  /6 

From  equation  (3.8.10)  it  follows  that 


(3.8.11)  / 


f(z) 


r*  (z-px) 


1 72" 


dz  =  - 


21 


4if(p1) 

i 


/pro 


/  y‘3/2[f(P1-P1y) 


-  f ( ) ]dy 


Consider 


(3.8.12) 


i  =  —  /  y'3/2[f(P,-P,y)  -  f ( Pi  )]dy 

0 


66. 


From  equation  (3.8.5)  we  have 


=  fCp^l-y)] 


=  a(n- 2)/2n  „,(n-2)/2  [1_pfJl(  1-y )]5/2C1+pP1(  1-y )]3/2 

Pi  ii-yj  « - ’17-5 — 

[l-P72(l-y)][p1(l-y)-p2]3/z 


and 


5/2, ,3/2 


f  Ox)  =  ^ 


.  _ ( n-2 )/2  *1_ppl*  (1+PPi> 


(1-p;;P1)(P1-P2)3/2  ’ 


so  that 


-  f(p1)  =  p[n'2)/2[(l-y)(n'2)/2g(y)  -  g(0)] 


and 


I  =  /1p[n_3)/2 y'3/2[(l-y)(n-2)/2g(y)  -  g(0)]dy  , 
0  1 


where 


[l-PMl-y)]5/2[l+pMl-y)]3/2 

g(y)  ~  5  t  3/2 

[l-p^l-y)]  [p1(l-y)-p2]'5/^ 


Expanding  g(y)  in  a  Maclaurin  series,  which  has  radius  of  convergence 
>  1,  we  obtain 


dy 


. 
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or 


(3.8.13) 


i  =  p["-3W[9(oivi§ii 

lx  “  i. 


where 


I o-l  y‘3/2[(l-y)(n‘2)/2  -  l]dy 


0 


and 


■  /V-3/2(i-y)(n-2)/2 

0 


dy  ,  for  k  =  1,2, 


Now  for  k  =  1,2, .. . 


1  ( k-1/2) -1  . 

I k  =  /  y  ( l-y)n/2_1dy  , 

K  0 


which  is  a  beta  function  with  parameters  (k-1/2)  and  n/2.  Thus 


(3.8.14)  Ik  ■  B(k-l/2,n/2)  •  %I^|nr)S(2) 


for  k  =  1,2,... 


We  may  integrate  by  parts  to  evaluate  Iq.  Letting 


u  =  (i_y)(n-2)/2  _  i  and  <j|v  =  y"3/2  ^ 


we  have 


f 


du  =  _  In^il  (i_y)(n-4)/2dy 


and  v  =  -2y 


-1/2 


so  that 


In  =  1 im 
u  £+0 


-2[(l-y)(n'2)/2-l] 


y 


VT 


1  .  (n-2)  /V1/2d-y)(n-4)/2dy 


0 


Now 


1  im 
£+0 


-2[(l-y)(n'2)/2-l] 

7^- 


i 

£ 


=  2  - 


1  im 
£+0 


-2[(l-e)(n"2)/2 

7* 


5 


which  upon  application  of  L'Hopital's  rule  becomes 


2  -  lim 
£->0 


(n-2) ( 1-e) (n"4)/2 


2  . 


Furthermore 


/  y-1/2(l-y)(n-4)/2dy  -  /y172-1!  1-y) ‘"'^^dy  , 
0  0 


which  is  a  beta  function  with  parameters  1/2  and  (n-2)/2.  Thus 


r1  -1/2,,  (n-4)/2  _  nr  i  ,o  ,woi  -  r(l/2)r[(n-2)/2] 

Jy  d-y)  °  B[  1/2 , ( n-2 )/2J  -  r[_i/2+(n-2)/2J—  ■ 


IQ  =  2  -  (n-2) 


r(l/2)r[(n-2)/2] 

rL(n-l)/2J  ' 


and 
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Recalling  that  r(a)  =  (a-l)r(a-l)  ,  it  follows  that 


(3.8.15) 


2  = 


on  r(l/2)r[(n-2)/2] 
l,W)  rL(n-l)/2J - 


2 ( n-2 )  r(l/2)r[(n-2)/2] 
2  rL(n-l)/2] 


o  r(l/2)r(n/2) 

L  r'LTn"-T)72J  * 


Using  equations  (3.8.14)  and  (3.8.15),  we  obtain 


T  _  r(l/2)r(n/2)  .  r(l/2)r(n/2)  -1  /T  ON 

A1  “  TLTn+r)/2'J  ITn-T)/2'Jrl'(n-l)/2j  "  TTTTT  1 l0"d)  ’ 


T  _  r(3/2)r(n/2)  _  ( l/2)r( l/2)r(n/2)  _  l  T  _  -l-l  /T  0n 

a2  "  r[Tn+3)/2] [(n+l)/2Jr[(n+l)/2J  "  Tn+TT  A1  "  rn-T)Tn+TT  li0"^ 


_  r( 5/2 ) r( n/2)  _  (3/2)r(3/2)r(n/2) 

3  "  TL(n+5)/2J  "  L(n+3)/2JrL(n+3)/2J 


3 

Tn+TT 


-1*1*3 

(n-1) (n+1) (n+3) 


n0-2) 


and,  in  general,  for  k  =  1,2,...  . 


T  _l*l*3***( 2k - 3 )  (  T  on 

Ak  "  ( n-l')'C n+TT(  n+3 )  •  •  •  ( n+2k- 3T  V  ^  ’ 


where  ( IQ-2) 


-2r( 1/2 )r( n/2)  _  -2/T  r(n/2) 
— rL(n-T)/'2'J  rL(n-l)/'2T  * 


From  equation  (3.8.13)  it 


follows  that 


70. 


(3.8.16)  I  =  p[n'3)/2{g(0)I0 

’  g<k>(0)  (-D  *1*3. . .( 2k-3)(  Iq-2)  ( 

+  k;1  ~T!  (n-l)  ( n+l )  Cn+3) . . .  (n+2k-3)  '  ' 

Combining  equations  (3.8.5),  (3.8.11),  (3.8.12),  and  (3.8.16)  we 
obtain 


h(  r) 


(l-p2)^  r4if(Pl) 


T 


+  2i 


4 up  /nr 
(n-3)/2 


{ 


P 


1 


'h 

'  »  nOOfm  ( -1) -1-3  —  (2k-3) ( In-2) 

9<0)  I„+  l  V1 

‘  u  k=l  K‘ 


(n-l)(n+l)(n+3) • ••( n+2k-3)/_ 


-  ^4~~ {-2pin-3)/2g(o> 


2itp  /nr 


1 


D 


+  B<n-3)/2|g(0)In  +  l 

u  k=l 


TV 


(-l)»l»3«*«(2k-3)(If)~2) 
( n- 1 ) ( n+l ) ( n-3 ) •••(n+2k-3) 


■  i(l-P2)1/2  p(n-3)/2[i  2]|g(0) 

2np3/nr  1  U  1 

”  g(k)(0)  (-1) >1-3 — (2k-3) 

kf1  F! -  (n-l)  ( n+l)  ( n+3)  - •  * ( n+2k-3") 


Since  g*k'(0)  depends  on  n  only  through  p,  ,  which  is  bounded 
between  0  and  1,  it  follows  that 


v  g(k)(0)  (-1) *1*3»..( 2k-3 )  .  n,„-l 

k|1  IT! -  (n-l) ( n+l)  (n+3)  •  ••(n+2k'-'jy  '  uk" 


■ 
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Furthermore,  in  approximating  |B|  we  incurred  errors  which  are 
relatively  CKn"1).  Thus 


where 


01/2 

h(r)  -  —  ~3PJ_ —  pjn‘3)/2[l0-2]g(0){l  +  0{ n “ 1 ) }  , 
2 up  /nr 


(i-PPl)5/2(i+PPl)3/2  (i-pPl)5/2(i+pPl)3/2 

9(0)  ■— 2;v,;  .\m  m —xwrzivr 


( 1“P  P  ^ )  ( ^l-^2^ 


i ( 1-p  f3  ^ )  ( ^2“^ i ) 


and 


pT  _oi  _  2/jF  r(n/2) 
Li0^J  "  ~  rL'(h-'l)72J  * 


Noting  that 


-3/2  _ 


1 


\-3/2  /I  -  P"f3 


-  P- 


202\  -3/2 


3„  3/ 2 


P  P 


1 


P2PX 


we  see  that 


g  ( o )  =  - 


p3p3/2(l-pg1)5/2(l+pP1)3/2  P3p3/2(1-Pp1) 

i(l-P2p1)(l-p2p2)3/2 


m-p^p 


Therefore,  for  0  <  Ipl  <  1,  we  have  the  following  first  approximation  to 
the  probability  density  of  r  : 


(3.8.17)  h(r )  - 


( 1-p2 ) 1/2  r(n/2) 


/nit 


rL(n-l)/2T 


pj72  (1-ppp 

- /—  {1  +  0(n  3)}  , 

_  /r  ( 1-p^p 1) 
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where 


 n( 1+p2 )  +  r(l-p)2  -  /fn 


(1+p2)  +  r(l-p)2]2  -  4n2p2 


2np‘ 


In  this  approximation  for  h(r),  each  omitted  term  involves  a 
factor  of  p  .  Also,  by  twice  applying  L'Hopital's  rule  it  may  be  shown 
that 


lim  p. 

p+0  1 


n 

(n+r)  * 


It  follows  that 


lim  h(r) 

p>0 


r(n/2) 

/ru  r[(n-l)/2] 


/r 


r(n/2) _ 1 

/mir  r[(n-l)/2]  ( 1+r/n )n/2 


which,  as  shown  in  section  3.7,  is  the  exact  distribution  of  h ( r )  in  the 
special  case  p  =  0  . 


-i  — 
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3-9  A  Renormalized  Asymptotic  Approximation  for  h (r).  The  Approximate 

Moments  of  h(r). 

From  the  approximation  (3.8.17)  we  see  that  for  0  <  |p|  <  1 

p?/2  (1-pP,)  , 

(3.9.1)  h(r)  =  Kn  -i- - J—  (1  +  CKn'1)}  , 

/F  (1-p2^) 

where  Kn  is  a  normalizing  constant.  The  value  of  Kn  may  be 
determined  by  studying  <j>(p^)  >  the  probability  density  of  . 

By  equation  (3.8.3)  we  have 


Pi  = 


n(l+p2)  +  r(l-p)2 


>J[n( 1+p2)  +  r(l-p)2]2  -  4n2P2 


2np‘ 


and 


n  n(l+p2)  +  r(l-p)2  +  ^[n(l+p2)  +  r(l-p)2]2  -  4n2p^ 

?2  ^ 


so  that 


n( 1+p2)  +  r(l-p)2 

Pi  +  P2 

np 


But  since  p2  =  (p2^)"1  ,  we  also  have 


Pi  +  P2 


P  P] 


+  1 


P  Px 


It  follows  that 


' 
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nP2(p1+p2)  -  n(l+p2) 

r  =  - - - 

(1-pT 

which,  upon  simplification,  becomes 


n(l+p2)]/l-P)2  , 


(3.9.2) 


r  - 


n(l-p1)(l-p‘;p1) 


(1-p)^! 


Furthermore, 


dr  C-nd-p2^)  +  n(l-g1)(-p2)](l-p)2p1  -  n(l-p1)(l-p2p1)(l-p)2 
^1  (1-P)4p2 


which  reduces  to 


(3.9.3) 


dr  -n ( 1-pB 1 ) ( 1+pB  x ) 
^1  (l-p)2|32 


The  probability  density  <t>( p is  related  to  h(r)  by 


(3.9.4) 


)  =  h(r ) 


dr 

n(  1-pp^)  (  1+pp^) 

h 

-nd-p^d-p^^)- 

dPi 

(i-p)2?f 

1 1 

1 — 1 

CEL 

CVI 

CL 

1 

rH 

_ 1 

for  0  <  p1  <  1  - 

Combining  equations  (3.9.1)  and  (3.9.4)  we  obtain 

-1/2 


n(l-pp1)(l+pp1)  ^rnd-^Ml-p  Pl^ 

*(Pl)=7T^  Vi 


< i-ppi  i  ,  ,i , 

- {1  +  0(n  1)} 

n-p\) 


* ,  1*  *  «"-‘u 


-n  {1.t 
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or 


(3.9.5)  *(p.)  = 


/n 


1'  "  TTTpT  Kn 


( l-pp1-p2p|+p3i3^)p[n-3)/2 
( l-p1> i/2( l-p2p1)3/2 


(1  +  0(n-1)} 


for  0  <  p^  <  1  . 
Consi der 


1  n  2  2  3  3%  (n-3)/2 

1  (1-pp^-p  (3^+p  p^)p^  2  3 

J  ~  J  "  i  /?  o  TT7)  “  P^p  _  P  ^3  +  P  J/i  > 

o  n-$l)iU(i-9dvl)6U  1  1  1  6  4 


where 


Ji  -  / 


(n-3)/2 

1 


o  (i-p1)1/2(i-P2p1)3/2 


dp. 


P 


( n-l)/2 
1 


/  - - rj7 - 2 -  ^1 

0  (i-p1)i/^(i-p^p1) 


1 


(n+l)/2 

J3 =  0J  (i-V/^T^V77  dP; 


,  and 


1  ,!"«>« 

J> '  I,  "1 ' 

From  equation  3  of  section  3.197  on  page  286  of  Gradshteyn  and  Ryzhik 
(1980),  we  have,  for  Re(a)  >  0,  Re(b)  >  0,  and  |d|  <  1,  that 


(3.9.6) 


ya-1(l_y)b_1(l_dy)"cdy  =  B(a,b)F(a,c;a+b;d) 


where  B(a,B)  is  the  beta  function  and  F(a,p;y;z)  is  the 


, 


' 


. 


hypergeometric  series 
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F(a,p;y;z)  = 


Thus 


Furthermore 


1  +  aPz  +  a(g+l)p(p+l)z2 
Y*  1  y( y+1 ) • 1*2 

(a+1)  (a+2)(3(  (3+1)  (  p+2)z3 

y(Y+l ) ( y+2 ) • 1*2*3 


n-1  l\.  r[(n-l)/2]r(l/2)  .  /Sr[(n-l)/2] 

2  ’  z)  fTn72T  fWD  ’ 


□  [n+1  l\  r[(n+l)/2]r(l/2)  .  (n-1)  /ir[(n-l)/2] 

B\T“  ’  2/  rL(n+2)/2J  n  r(n/2)  ’ 

0( n+3  l\  r[(n+3)/2]r(l/2)  .  (n-l)(n+l)  /fr[(n-l)/2] 

aYT~  ’  7)  rLn+4)/2J  "n(n+2)  r(n/2) 


and 


' 


n+5  1  _  r[(n+5)/2]r(l/2) 

2  ’  2  fL(n+6)/2] 


-  (n-1)  ( n+1)  ( n+3)  /iir[(n-l)/2] 
n( n+2) ( n+4)  r(n/2) 


From  equation  (3.9.5)  it  follows  that 


(3.9.7) 


UlPlj-1 

v'n 


=  (1-p)  r( n/2 )  J 3 
/n  /iir[(n-l)/2]  L  '2 

p(n-l)  r/3  n+1  n+2  2\ 

-  n  h\2  *  2  ;  2  ;  p  ) 

2  ( n-1 )  ( n+1)  Jz  n+3  n+4  2\ 

‘  p  n( n+2)  ^2  ’  2  ’  2  ’  p  / 


3  (n-1) ( n+1) ( n+3)  r/3  n+5 

p  n( n+2) ( n+4)  h \2  ’  2 


Note  that  for  p  =  0  , 


=  r( n/2 ) 

/nit  F[(n-l)/2] 


which  is  the  normalizing  constant  of  the  exact  distribution  of 
the  special  case  (p=0) 

Let 


h*(  r ) 


/r  (l-p2^) 


(3.9.8) 


78. 

It  follows  from  equation  (3.9.1)  that  the  moments  of  h*(r)  are 
asymptotic  approximations  for  the  moments  of  h(r).  Thus 


E(r)  »  /  rh*(r)  dr 
0 


1 

=  /'  rh*(r) 
0 


dr 

dp. 


which, by  equations  (3.9.2),  (3.9.3),  and  (3.9.8)  becomes 


E(r) 


1|"  nfl-^Xl-p^Pj) 

(1-p)2?! 


1/2 


n/2  (1  -pp^)  n( 1-pp^ ) ( 1+pp^ ) 

Pi  5  o  5  ^(3 

(1-P  PX)  (1“P)  Pi 


n3/2K  1  p[n-5)/2(l-p1)1/2(l-pg1)2(l+Pp1) 

- t/  - h - TTT1 - 

(1-p)3  0  (l-p^)  7 


n3/2K 


n 


(1-p) 


/ 

0 


p{n-5)/2(1-p1)l/2(l-pg1-p2Pi+p3p3) 


dp 


By  equation  (3.9.6)  we  obtain 


In  general  we  see  that 
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E ( rJ  )  «  /  rJh*(r)  dr 
0 


1  . 

/  rJh*(r) 
0 


dr 


dp. 


dp. 


=  K  / 
n  0 


■n(l-p1)(l-p2f!1)'|(2j'1)/2_n/2  ( 1-Pp1)  nd-p^Hl+pPj) 

-  5 - 5-75 - 


(1-p)  P 


1 


( 1-p231)  (1-P)2p2 


n 


(2j+l)/2  1 


(l-p)2j+1  0 


/  p[n-2j-3)/2  (1.Pi)(2J-l)/2  ( 1_p2p  , (2J-3J/2 


•  (l-pp1-p2p2+p3p'i1)dp1 

and  so,  for  j  satisfying  n  -  2j  -  3  >  0,  we  have 


(3.9.9)  E(r>) 


(2j+1)/2|< 

n 

tjtt 


(i-p) 


‘/n-2j-l  2j+l  V/-2J+3  n-2j-l  .  n  .  2^ 

LBV — Z—  ’  -T-)\-7-  ’  >2 '7 


-  p  B 1 


-2j+3 

n-2j+l 

n+2 

• 

v  2 

5  2 

5  2 

-2j+3 

n-2j+3 

n+4 

• 

5  2 

’  2 

/-2j+3 

.  n-2j+5 

n+6 

• 

3 • 10  The  Approximate  Probability  Density  and  Moments  of  F 

By  equation  (3.5.1)  we  see  that  the  statistic 


F  = 


n 


has  probability  density 


9(F)  =  T7TTJ"  h 


nF 


rfm 


From  equation  (3.9.1)  it  follows  that 


nK 


9(f)  =-^-e?/2 
( n-l)  1 


nF 

T7TTT 


-1/2  ( 1-pPi )  , 

- ^ -  (1  +  0(n_1) } 

(1-P  PP 


where  Kp  is  given  by  equation  (3.9.7)  and 


0  _  n(l+p^)  +  r(l-p)^  -  /[n(l+p^)  +  r(l-p)^]^  -  4n^p^ 

pi - 2 - 

1  2np^ 


n( 1+p^)  +  nF ( l-p)^/(n-l)  -  /fn(l+p^)  +  nF(l-p)^/(n-l)]^  -  4n^p^ 


Thus,  a  renormalized  asymptotic  approximation  to  the  probability  density 
of  F  is 


j— —  K„pT/2(1'pM  1 

(3.10.1)  g(F)  =  — — 5 — —  l1  +  °<n  >} 

’(n'1)  /F  (l-p2^) 


where 


■ 
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„  =  (1-p)  r(n/2) 

n 


/n7it  r[ ( n-1 ) /2]  l 


n-1 


n_ 

2 


p(n-l)  c/3 
n  ’ 

p2(n-l)(n+l)  J 
n(n+2)  F\ 


n+1  n+2  2  \ 

~~2~  *  ~7T  *  p 

3  n+3  .  n+4 

1  ’  “7“  »  IT  ;  p 


p3(n-l)(n+l)(n+3)  r/3  n+5  n+6 

n(n+2')'('n+4)  ’  T"  ’  ~T 


and 


Since 


e{fJ)  = 


( n-l) 


n 


E(rJ) 


we  obtain,  from  equation  (3.9.9),  that  the  jth  moment  of  F  is  given  by 
n1/2(n-l)jK  r 


E(FJ ) 


n 


( 1-p) 


Zj +T 


jn-2j+l  2j+l^  p/-2j+3  n-2j+l  .  n+2  .  2) 

-  pB\— 2 -  .  —£-)  Fy — 5-  .  - j—  >  ~T  ’  P  / 

2„/n-2j+3  2j+lV/-2j+3  n-2j+3  .  n+4  .  2^ 

-  P  B\  2  ■  t/v  ^  ^  ’  P  > 

,  3„/n-2j+5  2j+l\F/-2j+3  n-2j+5  .  n+6  2^1 

+  p  B^-j—  ,  —  jFlv  2  •  2  •  2  >  P  )}. 


for  j  such  that  n  -  2j  -  3  >  0  . 


■ 
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3.11  Areas  for  Further  Research 

By  analogy  to  Daniels'  (1956)  work  on  the  approximate  distribution 
of  serial  correlation  coefficients,  it  may  be  possible  to  refine 
approximation  (3.10.1)  so  that  the  error  is  relatively  0(n"3/2).  It 
may  also  be  possible  to  bound  the  terms  omitted  in  approximating 
g(F).  The  expression  for  Kn  could  be  put  in  a  series  form  and  studied 
further. 

For  selected  values  of  n  and  p  ,  numerical  integration 
techniques  could  be  used  to  obtain  the  approximate  distribution  function 
of  F.  This  could  then  be  compared  with  a  computer-simulated  empirical 
distribution,  and  with  the  limiting  distribution  of  F. 


* 
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APPENDIX  I 


EVALUATION  OF  A  KEY  DETERMINANT 
We  wish  to  evaluate  a  determinant  of  the  form 


(1.1) 


A  = 
n 


f  b 
b  a 
b 


(c) 


(c) 


a 

b 


b 

f 


The  following  method  of  evaluating  An  is  based  on  the  work  of  Dixon 
(1944)  and  Patton  (1961). 

Replacing  f  by  [a+(f-a)]  in  (1.1),  we  obtain 


A  = 
n 


a+( f-a) 
b 


b 

a 

b 


(c) 


(c) 


a 

b 


a+(f-a) 


=  (f-a) 


a 

b 


b 

a 

b 


(c) 


b 

a 


(c) 


a  b 

b  a+(f-a) 


n-1 


b 

a 

b 


(c) 


b 

a 


(c) 


a 

b 


b 

a+( f-1) 


n 
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Si  mi  1  arly , 


b 

a 

b 


(c) 


b 

a 


(c) 


a  b 

b  a+(f-a) 


n 


( f-a) 

a  b 
b  a  b 

(c) 

+ 

b  a  b 

•  • 

• 

(c) 

b  a  b 

b  a 

n-1 

b 

a  b 
b  a 


(c) 


(c) 


a 

b 


b 

a 


Thus, 


(1.2) 


An  =  Bn  +  2(f-a)Bp_1  +  (f-a)\_2 


where 


B  = 
n 


b 

a 

b 


(c) 


b 

a 


(c) 


a  b 
b  a 


Now, 


. 
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B  = 

1 

1 

1  . 

•  • 

1 

= 

1 

1  1  .  .  . 

1 

n 

0 

a 

b 

(c) 

-c 

a-c  b-c 

(0) 

0 

• 

b 

a 

• 

b 

•  • 

"C 

• 

b-c  a-c  b-c 

•  •  • 

• 

• 

(c) 

•  • 

b 

• 

a 

b 

• 

• 

(0)  b-c 

• 

a-c  b-c 

0 

b 

a 

n+1 

-c 

b-c  a-c 

■cR 

1+R2* 

-cR. 

L+R3»** 

•’"cRl+Rn+l 

Letting 

d  =  a 

-  c 

and 

e  = 

b  - 

c  ,  we 

obtai n 

Bn  = 

d 

e 

+ 

0 

1 

1  .  .  . 

1 

n 

e 

d 

e 

(0) 

-c 

d 

e  (0) 

e 

d 

• 

e 

•  • 

-c 

• 

e 

d  e 

•  •  • 

(0) 

• 

•  • 

e  d 

e 

• 

. 

(0) 

•  •  • 

e  d 

e 

' 

e 

d 

n 

-c 

e 

d 

n+1 

Hence 


(1.3) 


Bn  Bn  “  cDn+l 


where 


Cn  = 


e 

d 

e 


e 

d 


(0) 


n+1 


(0) 


e  d  e 
e  d 


n 


1 

and 
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D  = 
n 


0 

1 

1 


1 

d 

e 


(0) 


1 

e 

d 


(0) 


d 

e 


e 

d 


Now 

Cn 


=  d 


d  e 

-  e 

e 

e  d  e  (0) 

e  d  e  (0) 

e  d  e 

•  •  • 

e  d  e 

•  •  • 

(0)  e  d  e 

•  •  • 

(0)  e  d  e 

e  d 

n-1 

e  d 

n-1 


dCn-l  ‘  e  Cn-2  ’ 


or 


( I  -  4 )  Cn+2  ”  dCn+i  +  _  0  • 

Letting  v  be  the  forward  difference  operator  of  the  calculus  of 
finite  differences,  defined  by 

v[f(j)]  =  f(j+l)  , 


equation  (1.4)  may  be  expressed  as 
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(1.5)  (v2-dv+e2)C  =  0  . 

n 

The  auxiliary  equation  of  (1.5), 

x2  -  dx  +  e2  =  0 


has  roots 


u  = 


d+/?~ 


4e‘ 


and  v  = 


l-Zd2"- 


4e 


1 


Since  (1.5)  is  a  second  order  homogeneous  difference  equation  whose 
auxiliary  equation  has  distinct  roots,  u  and  v,  it  follows  that  the 
general  solution  of  (1.5)  is 


(1.6) 


cn  =  K1un  +  K2vn  . 


Using  the  initial  conditions 

Cx  =  K^u  +  K2v  =  d  ,  C2  =  Klu2  +  K2v2  =  d2  -  e2  , 


we  obtain 


,2  2  2 
d  -e  v 


2  2  2 

dv  -vd  +ve 

2  2 

uv  -u  v 


2  2 
dv-d+e 

u(v-u) 


and 
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d 


j 2  2 

d  -e 


.2  2  2 . 

ud  -ue  -u  d 

2  2 
uv  -u  v 


2  2 

d  -e  -ud 
v(v-u) 


Noting  that 


dv  -  + 


2  2 
e  =  -u^ 


and 


d 


5 


we  have  that 


d.7)  ■  -n^vr  and  k2  =  tc^vt  • 

Combining  (1.6)  and  (1.7)  we  obtain 

n+1  n+1 

(1*8)  ^n  ( u-v ) 

Now,  expanding  Dn  by  the  last  row,  we  obtain 


D„  =  (-1) 


n-1 


1 

e 

d 


(0) 


1  .  . 
d 


d 

e 


(0) 


e 

d 


n-1 


-  e 


0 

1 

1 


1 

d 

e 


1  (0) 


1 

e 

d 


d 

e 


(0) 


e 

e 


n-1 


+  d 


1 

d 

e 


(0) 


1 

e 

d 


(0) 


d 

e 


e 

d 


n-1 


or 


(1.9) 


Dn  =  (-1) 


n-1, 


■n-1 


-  eF 


n-1 


+  dD 


n-1 


where 


En  ’ 


1 

d 

e 

and  F  = 

0 

1 

1  . 

•  • 

1 

e 

d  e 

(0) 

n 

1 

d 

e 

(0) 

1 

e  d 

e 

1 

• 

e 

d 

e 

• 

• 

(0) 

•  •  • 

e  d 

• 

• 

(0) 

• 

e  d  e 
e  d 

1 

e 

n 

1 

e 
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Expanding  Fn  by  the  last  column,  we  get 


F  =  (-1) 


n-1 


n 


d 

e 


(0) 


(0) 


d 

e 


e 

d 

e 


n-1 


+  e 


1 

d 

e 


(0) 


1 

e 


(0) 


d  e 
e  d 


n-1 


or 


(I. 10) 


Fn  ■  f-DH’Vl  +  eD 


n-1  • 


Combining  equations  (1.9)  and  (1.10),  we  obtain 


Dn  =  (-1>n'lEn-l  *  et(-l)n-£En.rt-2]  +  dDn-l 


n-2, 


or 


(1.11) 


Dn  -  dDn_1  +  e^Dn.2  =  (-l)n-1(En_1+eEn_2) 


Expanding  En  by  the  last  row,  we  obtain 


. 


d 

e 
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En  -  (-D 


n-1 


e 

d 

e 


(0) 


(0) 


d 

e 


e 

d 


n-1 


+  e 


d 

e 


(0) 


e 

d  e 
e  d 


(0) 


d 

e 


n-1 


or 


(1.12) 


Ep  (-1)  +  eEn_i 


From  equations  (1.8)  and  (1.12),  we  get 


(1.13) 


^n  “  e^n-l  (“D 


,  n  n 
n-1  u  -v 

(u-v) 


Equation  (1.13)  may  be  expressed  as  a  second  order  difference  equation, 
the  solution  of  which  is 


(1.14) 


2e+d  + 


e(un-vn)  +  un+1 


(u-v)  (2e+d) 


From  equation  (I. 11)  it  follows  that 


°n+2  -  dDn+l  +  e\  “  <-Dn'1(WeEn)  , 


■ 
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or 


(1.15) 


(V2-dV+e2)Dn  =  (-l)n_1(En+1+eEn)  . 


From  equation  (1.14)  it  may  be  shown  that 


(1.16) 


r  .  „r  -  2en+1+(-l)n(un+1+vn+1) 
n+1  eLn  Ze+cl 


Combining  equations  (1.15)  and  (1.16),  we  obtain 


(1.17) 


oooi  a\n+l  /..n+ln+L 

[^Wl»„  ■  a-«l  ;<»  »■  ' 


To  obtain  the  complementary  function  of  equation  (1.17),  we  consider  the 
homogeneous  difference  equation 


(V2-dV+e2)Dn  =  0  . 


The  corresponding  auxiliary  equation. 


-  dx  +  e^  =  0 


has  roots  u  and  v.  Thus 


Dn  -  x  un  + 
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is  the  complementary  function  of  the  complete  solution  of  equation 

(1.17).  To  obtain  a  particular  solution,  we  assume  a  solution  of  the 
form 

(1-18)  Dn  =  P(-e)n+1  +  Qnun  +  Rnvn  . 

Then 

(v2-dv+e2)Dn  =  P(-e)n+3  +  Q(n+2)un+2  +  R(n+2)vn+2 

-  d[P(-e)n+2  +  Q(n+l)un+1  +  R(n+l)vn+1] 
+  e2[P(-e)n+3  +  Qnun  +  Rnvn]  , 

from  which  it  may  be  shown  that 

(1.19)  (v2-dv+e2)Dn  =  -P( 2e+d) ( -e)n+2  +  Q(u-v)un+1  -  R(u-v)vn+1  . 

Equating  coefficients  in  equations  (1.17)  and  (1.19),  we  obtain 

2-1  1 
P  ~  e( 2e+d) 2  *  ^  ”  (u-v)(2e+d)  ’  R  "  (u-v)(2e+d)  * 


and  equation  (1.18)  becomes 


96. 


Dn  = 


.  2(-e) 


n+1 


e( 2e+d) 


nu11  nv11 

(u-v) ( 2e+d)  ( u-v ) ( 2e+d) 


2(-l)n-1en 

(2e+d)2 


-  n 


un-vn 


(2e+d) (u-v) 


Thus  the  complete  solution  of  equation  (1.17)  is 

n  n  "I 
u  -v 

(2e+d)(u-v) 

where  we  may  use  the  initial  conditions  to  solve  for  x^  and  x2  • 
defi ni tion 


(1.20) 


D  =  X,un  +  X«vn  +  2(~1) 
n  1  2  ( 2e+d) 2 


ri-i  n 

e 


-  n 


D 


1 


1 0 1  =0  and 


5 


and,  from  equation  (1.20), 


x^u  +  \2V 


2e  (u-v) 

,2e+d)lu-vl 


X^U  +  X2V 


(u+v) 

( 2e+d)2 


and 


2  2 
D2  —  x^u  X2 V 


2(u2-v2)  .  ,  „2  t  ,  „2  _  (u2+v2)  _  , 

^7  -  ( 2e+d )  ( u^vT  ‘  V  V  ^7 


2e 


By  inspection  we  have 


h  - 


1 


( 2e+d) 


1 


and 


x2 


1 


(2e+d) 


so  that  equation  (1.20)  becomes 


(1.21) 


n  .  un+vn+2( -l)n’1en 
n= 


n( un-vn ) 
(u-v) ( 2e+d) 


Combining  equations  (1.3),  (1.8),  and  (1.21)  we  obtain 


un+1-vn+1  (n+l)  (un+1-vn+1) ) 

u-v  '  c\  (2e+d)2  -  (u-v ) ( 2e+d)  t* 


or 


(1.22) 


Bn  = 


u 


n+l 


(u-v)(2e+d) 


-j(2e+d)2  -  c[(u-v)-(n+l)  (2e+d)]| 


-v 


n+l 


( u-v) ( 2e+d) 


r<J(2e+d)2  +  c  [( u-v )+( n+l )  (2e+d)]j> 


2c (-e) 


n+l 


( 2e+d) 


Substituting  this  value  of  Bp  in  equation  (1.2)  we  find  that 


An 


u 


n+l 


(u-v)(2e+d) 


^-|(  2e+d) 2  -  cC(u-v)  -  (n+l)  (2e+d)]j> 


+  ^  f~alu -  i(  2e+d) 2  -  c[(u-v)  -  n(2e+d)]l 

( u-v)  (  2e+dr  '■ 


(f-a)  u 


2  n-1 


( u-v) ( 2e+d) 


n+l 


( u-v)  ( 2e+d) 


-j(2e+d)2  -  c[(u-v)  -  (n-1) ( 2e+d) ] j 
-|(2e+d)2  +  c[(u-v)  +  (n+l) (2e+d)]| 


■ 
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2( f-a)  v 


n 


(u-v)(2e+d) 


■^■|(2e+d)2  +  c[(u-v)  +  n(2e+d)]J> 


(f-a)  v 


2  n-1 


( u- v ) ( 2e+d ) 2  ^ 


i(2e+d)2  +  c[(u-v)  +  (n-1)  ( 2e+d) ij1 


.  2c( -e)n+1  A  4(f-a)c(-e)n  ^  2( f-a)2c( -e)n_1 

- 2“  + - 7777772  + - 777772 - 


( 2e+d) 


( 2e+d) 


( 2e+d) 


Collecting  terms  and  simplifying  we  obtain 


(1.23) 


An  " 


u 


n-1 


(u-v)(2e+d) 


7|[u+(f-a)]z[(2e+d)^  -  c{(u-v)-n(2e+d)}] 


+  c(2e+d)[u2-(f-a)2]} 


n-1 


(u-v) ( 2e+d) 


2  ^ [v+(f-a)]2[(2e+d)2  +  c{ ( u-v)+n( 2e+d) }] 


c( 2e+d)[v2-( f-a) 2] | 


+ 


o  /  \  n- 1 

2c( -e) 

( 2e+d ) 2 


Ce-( f-a ) ] 2 


where  u  and  v  are  roots  of  the  equation 


-  dx  +  e^  =  0 


and 


99. 


d  =  a  -  c  ,  e  =  b  -  c  . 


It  is  easily  shown  that  An  is  symmetric  in 


u  and  v. 


APPENDIX  II 


THE  JOUKOWSKI  TRANSFORMATION 
Consider  the  complex-valued  mapping 
(II. 1)  w  =  z  +  i  , 

i  0 

which  is  known  as  the  Joukowski  transformation.  Letting  z  =  re  we 
have  that 

_  i  Q  _  i  —  i  0  — 1 

w  =  re  +  r  e  =  r(cos  0  +  i si n  0)  +  r  ( cos ( -0 )  +  isin(-e)) 


or 

w  =  (r+r"^)cos  0  +  i  (r-r’^sin  0  , 
which  is  a  more  convenient  expression  to  study. 

If  r  =  1  then  w  =  2cos  0  .  Thus,  as  z  proceeds 
counterclockwise  along  the  unit  circle  from  (1,0)  to  (-1,0)  to  (1,0), 
w  moves  along  the  real  axis  in  the  w-plane  from  2  to  -2  and  back  to  2. 

If  z  moves  on  a  circle  of  radius  r  ^  1  centred  at  the  origin, 
then  the  locus  of  w  =  u  +  iv  is  the  ellipse 


2  2 
u  +  v  =  1 
a  b 

with  a  =  r  +  r“^  and  b  =  r  -  r~^.  It  follows  that  if  |z|  =  r  >  1, 
then:  (i)  as  r*  +®  ,  a  and  b  -►  +»  ;  (ii)  as  r  ->  1,  a  -*  2  and 

b  -*  0;  (iii)  a  >  2  and  b  >  0;  (iv)  r  >  r“ ^  ;  (v)  for  0 

<  0  <  %  ,  w  lies  in  the  upper  half  of  the  w-plane;  and  (vi)  for 
n  <  0  <  2n  ,  w  lies  in  the  lower  half-plane.  Thus  the  Joukowski 
transformation  maps  the  w-plane,  cut  along  the  real  axis  from  -2  to  2, 
inclusive,  onto  the  exterior  of  the  unit  circle  |z|  =  1.  Furthermore, 
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■ 
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the  upper  half  of  the  w-plane  is  mapped  onto  the  upper  half  of  the 
exterior  of  l z I  =  1,  and  the  lower  half  onto  the  lower  half  exterior. 

Similarly,  if  |z|  =  r  <  1,  then:  (i)  as  r  -►  0,  a  -►  +=°  and 
b  -►  -®  ;  (ii)  as  r  +  1,  a  ■>  2  and  b  >  0;  (iii)  a  >  2  and  b  <  0; 

( i v )  r  <  r“*;  (v)  for  0  <  9  <  %  ,  w  1 ies  in  the  lower  half  of  the 
w-plane;  and  (vi)  for  n  <  0  <  2n  ,  w  lies  in  the  upper  half-plane. 
Thus,  the  Joukowski  transformation  also  maps  the  cut  w-plane  onto  the 
interior  of  the  unit  circle  |z|  =  1,  with  the  upper  half  of  the  w-plane 
going  to  the  lower  half  interior  and  the  lower  half  of  the  w-plane 
mapping  to  the  upper  half  interior. 

We  may  also  express  equation  ( 1 1 . 1 )  as  the  quadratic  equation 

z^  -  wz  +  1  =  0 

with  roots 


Observing  that  z^z2  =  1  we  obtain 

I z1 | I z2l  =  1. 

It  follows  that  for  given  w  either  one  solution  of  equation  ( 1 1 . 1 ) 
lies  inside  the  unit  circle  and  the  other  outside,  or  both  solutions  lie 


on  the  unit  circle. 


■ 
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